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ABSTRACT 

The non-linear elasto-plastic response of a clamped 
flat plate to a typical air shock wave were investigated. 

The nonlinear effects to the plate responses due to the 
material and geometric nonlinearity were studied. 

In this study, (1) the necessity of the modification of 
old armored vehicles was reviewed and (2) the NASTRAN code 
was employed in this investigation, (3) the theoretical 
background of the nonlinear transient analysis was described, 

(4) a step by step procedure of analyzing the dynamic load 
problem by shock wave using NASTRAN code was developed, and 

(5) sensitivity analyses were performed and also difficulties 
associated with the nonlinear analysis were described. 
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INTRODUCTION 



I . 

Conventional tanks (old tanks) which have steel armor, 
can be completely defeated by one round of modern projectiles 
such as anti-tank missiles. That has been proved in the 1973 
Arab-Israeli war. In that war, the 40% of total Israeli 
casualties were caused by armored vehicle destructions. Thus, 
tank armor improvement became the first priority in the Israeli 
tank corps . 

If war occurs in any country, the same problem will be 
faced. That is the reason why the modification of the conven- 
tional tank armor is needed against modern anti-tank weapons . 

To improve the armor protection, the material properties, 
projectile characteristics, and explosion phenomena must be 
understood well. This is a complicated and difficult problem 
and is originally a nonlinear problem. 

First the shock wave from an explosion imposes a dynamic 
load on any object. That dynamic load is characterized by 
rapidly reached peaik overpressure which then decreases as the 
shock wave decays. The net effects of the load depend both 
on the nature of the shock wave and the geometry and construc- 
tion of the object. For example, if the explosion occurs on 
thin skirting plate, perforation takes place by impact load 
and their impulse loading for main armor would be reduced 
significantly. The typical air shock wave was characterized 
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as follows; the shock wave is supersonic and generates the 
abrupt increasing overpressure. 

The non-linear elasto-plastic responses of a clamped flat 
plate to a typical air shock wave were investigated. The 
nonlinear effects to the plate responses due to the material 
and geometric nonlinearity were studied. 

In this study, (1) the necessity of the modification of 
old armored vehicles was reviewed, (2) the NASTRAN code was 
employed, (3) the theoretical background of the nonlinear 
transient analysis was described, (4) a step by step proce- 
dure of analyzing the dynamic load problem by shock wave 
using NASTRAN code was developed, and (5) sensitivity analy- 
ses were performed and also difficulties associated with the 
nonlinear analysis were described. 
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II . LITERATURE REVIEW 



A. MODERN PROJECTILES 

An important version of the AP (Armor-Piercing) round is 
that known as APDS (Armour Piercing Discarding Sabot) . In 
this round the armour-piercing element is a hard core of 
calibre significantly smaller than that of the complete 
projectile and shaped for maximum anti-armour effectiveness. 
This is surrounded by a much lighter casing, or sabot, which 
breaks up as the projectile leaves the gun barrel. The 
point of this arrangement is that the core can have a shape 
which differs considerably from that of the complete projec- 
tile and, not having been scored by the rifling in the barrel, 
is aerodynamically clean. Moreover, since it is very much 
heavier than the sabot, the core acquire a greater energy 
than could be imparted to it if it were fired, without sabot, 
by a gun of appropriate calibre and the same chamber pressure 
as that of the larger gun. A controversial refinement 
involves the use of depleted uranium in the core. This 
metal has an exothermic reaction with steel which aids 
penetration; it is controversial because of its connection 
with the other isotopes of uranium and their radioactive 
properties . 

High explosive (HE) rounds with contact fuzes can also 
be used against tanks and may, if they are large enough, be 
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effective against exposed portions of the tank mechanism 
(such as the tracks) but will not in general be effective 
against heavy armour. Two modern refinements of the HE 
round may be effectively used in some circumstances: they 

are the HEAP (High-Explosive Armour Piercing) and HESH 
(High-Explosive Squash-Head) rounds. 

The HEAP round has a projectile of relatively high 
velocity, a hardened case and a high-explosive charge which 
is contact fuzed with a very short delay. The effect of the 
arrangement is that the projectile either embeds itself in 
the armour or passes through it and then explodes with an 
effect which is more damaging than that of the simple HE 
round . 

HESH rounds can be used with lower-velocity weapons, and 
are so designed that the nose of the shell crumples on impact 
to make a large contact area before the explosive is detonated 
Such rounds will not, in general, pierce armour, but the 
explosion sets up shock waves in the armour which in them- 
selves can incapacitate the crew and may also cause pieces of 
metal to break away from the surface remote from the explosion 
Such fragments fly off at high velocity and can do great 
damage inside a tank turret. 

Most widely used in infantry weapons, however, is the 
projectile (which may be a shell, a rocket, a grenade or a 
guided missile) with a shaped-charge warhead frequently 
described as a HEAT (High Explosive Anti Tank) round. This 
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type of projectile is designed to be effective with only a 
moderate impact velocity and depends for its effect on the 
design of the warhead and is formed round a hollow cone with 
its apex to the rear. The base of the cone is of approxi- 
mately the same diameter as the warhead at that point and in 
front of it is an air space enclosed by the rounded or 
pointed head of the projectile on which the fuze igniter is 
mounted. The operation of the fuze is critical to the success 
of the charge. It has to have a set distance from the surface 
of the armour when it is detonated and for this reason the 
igniter is carried well in front, often giving rise to the 
distinctive pointed nose of a hollow charge round. 





PARTICLES 



STEEL ARMODB 



Figure 2.1. The Effects of Current Anti-Tank Projectiles 

on the Steel. 



The igniter is connected to the detonator which is at the 
back of the charge, this starts a detonation wave running 
through the conical explosive and the result is that the 
force of the explosion is concentrated into a thin jet 
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moving at a speed of over 20^000 m/s and at a great tempe- 
rature. The jet takes a finite distance to form, hence the 
importance of having the igniter in front and the distance 
required to form a jet is known as the "stand off". For most 
modern warheads the stand off distance needs to be between 2 
and 3 diameters of the warhead and to provide the full dis- 
tance sometimes interferes with the ballistics of the 
projectile . 

The effects of a hollow charge on the interior of a 
vehicle can be variable. A very small jet may well not 
cause much damage apart from making a hole in the armour and 
destroying whatever is in its direct path, but it is possible 
to "turn" warheads so that the penetrative effect is lessened 
but there is more energy left after passing through the 
armour and this energy can be spread out into more of a fan- 
shaped exit blast with appropriately enhanced damage. All 
hollow charge jets carry a certain amount of molten metal 
into the interior and some research has gone into finding 
ways of carrying more through the hole. 

Two examples of the Soviet anti-tank guided missiles are 
described below. [Ref. 1] 

1. Sagger Anti-Tank Guided Missile 

Sagger is the NATO code-name for the Soviet anti- 
tank missile which is also known by the US alphanumeric 
designation AT-3, which is now extensively deployed in 
Europe and elsewhere. Currently it is known to be used in 
at least three ways: 
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a) As a man-portable missile 

b) Mounted on a bracket 



c) Mounted on a vehicle 

The characteristics of sagger missile is as follow. 



a) 


Type: Wire-guided MCLOS (Manual Command to Line 

of Sight) 


b) 


Weight at launch: 11.3 kg 


c) 


Diameter: 120 mm 


d) 


Warhead: HEAT 


e) 


Range: 500-3000 m 


f) 


Penetration: 400 mm 


2 . AT- 


-4 Anti-Tank Guided Missile 



There have been reports for some years of a Soviet 
SACLOS (Semi-Automatic Command to Line of Sight) system, but 
it was not until 1980 that photographs and a little informa- 
tion appeared in some Warsaw Pact military journals. 
Characteristics : 



a) 


Type: Wire guided SACLOS system 


b) 


Weight at launch: 40.4 kg 


c) 


Diameter: 120 mm 


d) 


Warhead: HEAT 


e) 


Range: 2000 - 2500 m 


f) 


Penetration: 500-600 mm 
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Figure 2.2. Sectioned View of Sagger Missile. 

Which of these various types of projectile will be most 
effective against armour which depends on the design and 
thickness of the armour. Against solid armour the HEAT 
round is particularly effective; but in certain circumstances 
a sandwich construction is likely to be more vulnerable to an 
APDS round. There is therefore a good deal to be said for 
having a mix of projectiles available for satisfactory 
defense against tanks. 

B. CONFIGURATION OF ARMORED VEHICLES 
1 . Armor 

We can define the military armor as a material 
intended to protect men and installations of equipment. 
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Aluminum, titanium and magnesium alloy, plastic laminates, 
nylon fabrics and masticstone composites have all been used, 
but the most common form of armor is heat-treated alloy steel. 
This comes mainly in two forms: homogeneous and face-hardened. 

The first has constant ballistic and mechanical quali- 
ties throughout, which in modern tank armor means high tensile 
strength (of the order of 150,000 lb sq in) and medium hard- 
ness (i.e. Brinnell N 300); face-hardened armor consists of 
a very hard first layer which resists penetration but lacks 
toughness and a more elastic core to support the face. Naval 
deck armor and modern tank armor are homogeneous, while face- 
hardened armor is still used for some light armored vehicles 
and side protection for naval vessels. 

Steel armour usually means rolled or wrought plates, 
but castings are used for tank turrets and even for entire 
tank hulls, as in the M.48 and M.60 tanks (Steel armor is 
usually a nickel-chrome-molybdenum alloy, but other alloys 
including manganese-silicon-chromium-nickel steel have been 
used, especially in Soviet tank armor.) . 

Typically, a well-designed hollow-charge warhead can 
be expected to penetrate armor up to five times the diameter 
of the warhead. The best defense against such attack is not 
to make thicker armor, but to dissipate the jet before it 
reaches the surface. To do this a light screen is carried 
some distance in front of the main armor, of sufficient 
thickness to set off the fuze of any hollow-charge projectile. 
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Quite naturally such screens do not last long, but they are 
easily replaced and provide almost complete protection for 
the armor behind. A common use of such a technique is the 
thin skirting plates carried over the tracks of all modern 
tanks. Before leaving this subject it is appropriate to 
mention a recent British armor development, the 'Chobham' 




Figure 2.3. The Effect of Anti-Tank Projectiles on "Cobham" 

Armo r . 

armor. This is claimed to offer greatly enhanced resistance 
to all forms of projectile as is indicated by the accompanying 
illustrations. Details of the construction of this armor have 
not been released, but it seems probable that it is a sand- 
wich made up of two sheets of armor plate separated by a 
highly dispersive, possibly granular, filling. The disad- 
vantage of this armor lies in its bulk. To gain the protec- 
tion, it has to be thick and fairly heavy; tanks are already 
large enough and the addition of armor plating, which is 200 
mm or more in thickness makes a startling difference to the 
outline and bulk of the entire vehicle. 
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2 . Current Modification of Old Tanks 

Determining which one should we choose, a new one 
or modification of an old one usually depends on cost and 
political situations. 

The U.S. completed conversion of M48 tanks to M48A5 
standards in 1980. This conversion basically included adding 
diesel engines to M48Als and adding 105-mm guns to the tur- 
rets of M48Als and M48A3s. Numerous minor changes were also 
made. But little improvement in armor protection has been 
made. If we want to take more modifications for M48A3s and 
M48A5s tanks, a spaced applique armor should be added and 
include a 1,200 hp engine, improved transmission, better 
fire control, and a hydropneumatic suspension system. 

The following example shows modern modifications. 

The High Performance M60 MBT has been developed as a 
private venture by the General Products Division of Teledyne 
Continental Motors. It is essentially an M60 series MBT 
fitted with a new powerpack, new suspension and, as an option, 
additional armor. These modifications can be carried out on 
existing M60 series. Teledyne believes that to fit an M60A1 
with the new powerpack would cost under $200,000, while the 
added armor would cost just $40,000. 
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SPACED APPLIQUE ARMOR 
EQUALS 40.64 CM RHA 




Figure 2.4. Applique Armor in M60 MBT . 

3 . Israel Add-On Armor 

Israeli experience in the 1967 campaign proved that 
mobility was not substitute for armor protection and they 
therefore decided at an early stage that the main emphasis 
would be placed on armor, with firepower and mobility second 
and third priorities. After that they spent a lot of money 
to figure out this protection problem. Total cost of 
research, development, trials and the building of the proto- 
type vehicles was about $65 million. 

In May 1977 Israel finally announced that it had 
developed a new MBT called the Merkava which was made of 
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cast and welded armor with a well-shaped glacis plate. 

Behind the first layer of cast armor is a space filled 
with diesel fuel and then another layer of armor. This 
spaced armor gives the tank protection from HEAT projec- 
tiles and ATGWs (anti-tank guided weapons) . 

They also modified their old tanks such as M60s, 
Centurions, etc. 

The add-on and readily detachable protective arrays 
seen on Israeli M60s and Centurions contain ceramic tiles, 
would also mark a first in terms of operational deployment. 
The possibility cannot be excluded that armor of both kinds 
is involved. 

While no answers are to be found in external appear- 
ance, the configurations, as seen in the accompanying Defence 
Attache photographs, show obvious differences in the thick- 
ness and shape of the packs provided for M60 and Centurion 
respectively. The Centurion's additional armor is mainly in 
the form of quite shallow panels fitted to the brow of the 
turret and thicker pieces to the glacis and the mantlet 
either side of the main gun; a hole is included in the left- 
hand piece to permit the firing of the coaxial machine gun. 
Wedge-shaped packs are mounted on top of the storage boxes 
fitted forward above the track guards. 
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Figure 2.5. 



Add on Armor Blocks 



in M60. 
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Figure 2.6. Add on Armor Blocks in Centurion. 

The M60/ partly because it lacks whatever protective 
values come with the Centurion’s storage bins, has been 
provided with a considerable amount of add-on armor around 
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the turret. Here the packs are thick and vertically deep, 
with smaller and thinner sections patched in — with no 
great concern for aesthetics — to the upper surfaces which, 
because of the turret's curvature and varying angles to the 
perpendicular, not to mention obstructions such as lifting 
eyes, cannot be covered by the main boxes. As with the 
Centurion, panels are fitted to the front of the turret roof 
and to the glacis; in the latter case, because the M60's 
basic glacis is curved, it has been necessary to fit a 
straightening frame to which the packs can be attached. The 
M60's installation is completed by the fitting of wedge- 
shaped or angled sections either side of the driver's hatch 
to help protect the vulnerable area where turret meets hull, 
with further angled sections each side of, and thinner plates 
above, the mantlet. 

Although it can be stated with a fair degree of 
certainty that the add-on arrays are intended to help defeat 
shaped-charge high-explosive anti-tank (HEAT) attack, any 
discussion of what the packs and panels consist of must, in 
the absence of any word of guidance from the Israelis - who 
on the contrary are happily watching, if not promoting, the 
circulation of a number of often fanciful and frequently 
mutually exclusive 'solutions' to the puzzle-fall short of 
a definitive answer. 

While one source suggested that the boxes were empty 
(and the large ones on the M60 do sound quite hollow when 
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thumped) the probability is that there is more to it than 
that. In view of solutions that have been adopted else- 
where, though usually in the context of new-build tanks, it 
seems quite likely that the Israelis have developed their own 
ceramic tiles and that the boxes are partially filled with 
them. The alternative of a form of active armor is suggested 
more by the vehemence with which several sources in Israel 
stated this to be the case than by any objective evidence. 

A number of different types and configurations of 
active armor are possible. An example is a sheet of explosive- 
laid of top of a steel plate, with a covering for environmental 
protection on top - which detonates when hit by the incipient 
jet of a shaped-charge warhead and in the process disrupts the 
jet's further formation. If the Israelis are indeed deploying 
a form of active armor, such a fact would carry the implica- 
tion that they have sufficiently overcome the formidable 
difficulties that seem to have prevented larger and better- 
endowed countries from doing likewise; the United States and 
Soviet Union are among the nations that are known to have been 
engaged in related research and experimentation for many years. 

C. ELEMENTARY DETONATION THEORY 

1. The following assumption will be taken for the 
simplest analysis of a detonation: [Ref. 2 ] 
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a) The flow (motion of the shock front) is one 
dimensional 

b) The plane detonation front is a jump discontinuity 
at which chemical reaction is assumed to be 
instantaneous . 

c) The discontinuity is time-independent, so that the 
product composition (at least directly behind the 
shock front) is likewise time-independent. 

u represents the particle velocity of blast wind, which is 
that of the products behind the shock front. Subscript x 
represents the initial values of the properties, y those 
after shock. D is the detonation speed and p is the shock 
wave overpressure. 

From mass conservation across the front we can write 

P D = p (D-u ) (1) 

X y p 




Figure 2.7. Explosive Shock Front. 
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and from conservation of momentum, 

p-p=puD (2) 

X p ' 

Eliminating up from the two equations we obtain the equation 
for the Rayleigh line; 



- (py - p^)/(v^ - Vy) = 0 



(3) 



(Observe that in part of the previous equation we replace 1/p 
by the specific volume, v) . Since for a detonation, p > p / 

y ^ 

it follows that v < v (for real values of the detonation 

y X 

velocity) ; and, indeed, it is observed as a rule that there 
is about a one-third increase in density on detonation. If 
this condition is substituted in eq (1) , it follows that 
0 < Up < d, so that the products travel toward the detonation 
front, though trailing behind. 

If the Rayleigh equation is solved for p^ we obtain 



p = (p + P D^) - Px^D^ V 

^y ^x X y 



(4) 



wherein p^ is linear with v^, the slope of the line being 
2 2 

-P D . The conservation conditions thus require that the 
initial and final states of a detonation must lie at the 
termini of the Rayleigh line. 

Conservation of energy leads to 

E(P^ , X=l)+P v^+l/2(D-U = E (P^,l/x,X=0)+P V +1/2D^ (5) 

yy yy p 

in which X = the extent of chemical reaction (X =0 for no 
reaction and 1 for complete reaction) . Eliminating u and 
D from Eq (5) , suing Eqs (1) and (2) leads to the Hugoniot 
curve : 
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E(P ,v ,A=1) 

y y 



( 6 ) 



- ‘VV <V'^y>=° 

A Eh = 1/2 (Py+P^) (Vx-Vy) (7) 

AEy is the internal energy of the detonation products at the 
very high pressure and temperature of the detonation state 
less the internal energy of the unreacted explosive at ambient 
temperature and pressure. It may be written as the sum of two 
terms : 



AEh = AE^ + AEp 



( 8 ) 



where AE^ is the internal energy change for the conversion 
of reactants to products isothermally , and AE^ is that 
required to warm and compress the products from ambient p 
and T to the detonation values. AE^j, is often symbolized 
as -Q in explosive literature (Q being the heat evolved in an 
explosion when corrected to isothermal conditions) . If the 
products were ideal gases, then AE^ would be given /c^dT, 
and this approximation is often made, though the products at 
pressures of hundreds of bars, are far from being ideal gases. 

2. The Hugoniot curve is that of a hyperbola in the p-v 
plane. For the non-reacting case (A=0) the hyperbola passes 
through the initial state, p ,v ; but for the case we are 
considering here (A=l) it does not. The conservation laws 
require that the final state, however, lies on the Hugoniot 
curve as well as on the Rayleigh line. We consider below 
several cases: 
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For a too small D(=d 2 ) the Rayleigh line does not 
intersect the Hugoniot and thus no solution exists. For 




Figure 2.8. The Hugoniot Curve. 

large D (=Dj^) , there are two intersections, called strong 

(labelled S) and weak (W) . For one particular value of 

D (=D .) the Rayleigh line is tangent to the Hugoniot. 

^ 1 

It can be shown that the relative particle velocity, 

D-u , is subsonic at S (D=u < a ) , supersonic at W and 
P P ^ 

sonic at the point of tangency, called the Chapman-Jouguet 
point: (D-Up)^j = a^ (sonic velocity in the shocked medium) . 

There exist a number of computer codes based on such 
non-ideal gas treatments (RUBY and TIGER are two of these) and 
moderate success has been achieved in computing detonation 
properties from these codes. The equations are; 
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(9) 



Py = kP^2 ^ 

D = A4)^ (1 + BP^) (10) 

4) = NM^Q^ (11) 

with k = 0.762, A = 22.3, B = 0.0013 (SI units) . Q, the 
heat of explosion in J/kg; N, the number of moles of gas per 
kg; and M, the average molar mass of the gas (kg/mole are all 
found from an assumed stoichiometry of decomposition, combined 
with thermochemical principles. Although N, M, and Q are all 
strongly affected by the assumption made as to stoichiometry, 
it turns out that the important parameter, 4> , is rather 
insensitive to the assumption. 
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III. ANALYSIS OF DYNAMIC RESPONSE OF PLATE BY SHOCK WAVES 



A. THEORETICAL BACKGROUND 

1. Nonlinear Static Analysis 

The fundamental method used to obtain a static solu- 
tion is to minimize the error vector {6}, given by the 
following equation: [Ref. 7J 

{6} - (P}+ {Q}-{F} (12) 

where {6} is the error vector of unbalanced forces acting 
at all grid points, (P) is the vector of applied external 
loads, which may change with displacements, {Q} is the unknown 
vector of constraint forces due to single and multipoint con- 
straints, and {f} is the vector of grid point forces generated 
by element motion and stress. The terms are functions of 
displacement, temperature, and stress history. 

It must be noted that degrees of freedom not involved 
in constraints produce null terms in {Q} and that dependent 
constraint points produce no errors. Therefore, when Eq (12) 
is reduced to the solution coordinates the constraint forces, 
{Q}, disappear. Terms in the vector {F} are dependent on the 
deformations of the finite elements. In linear static analysis 
the vector {F} becomes 

{F^inear> = [K^lfu} - {p'^(T) } (13) 
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where [k^] = The linear stiffness matrix 



T 

{P } = The "thermal load" vector, i.e., a 

vector of grid point which produces the 
same element strains as the temperature 
distribution , 

{u} = The vector of grid point displacements. 

To obtain a nonlinear solution, a "tangent" matrix 
[K] is used in which the terms are derivatives in the form 



K 



6Fi 

ij ^ ^6U. = u 



(14) 



The matrix [k] contains both geometric and material non- 
linear effects. The approximation occurs because the forces 
are highly variable with displacements and nonsymmetric terms 
are approximated. If {u^} is a known displacement which 
produces an error {6(u^)}, then Newton's method may be used 
to predict a new vector with a smaller error. For 

a displacement {u} near a known solution, {u^}, the non- 
linear force is approximately 

{F(u)} (F(u^)} + fK]{u-u^} (15) 

The substitution of Eqn (15) into Eqn (12) with {6} = 0 
yields the result 

[K]{u-u^} {P} + {Q} - {F(u^)} (16) 

which can be written in terms of the error {6} to provide 
the format of the Newton-Raphson iteration method: 



[K] {u^‘*’^-u^} = {6 (u^) } 



(17) 
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When [K] may be inverted, Eqn (17) may be used to solve 
for = {(u)}, a new estimate of the solution. 

Equation (17) provides an incremental form for solution 
iteration. Another form of the iteration equation is 
obtained from Eqns (12) and (17) ; 

[K] {u^‘‘'^-u^} = {p^ + p"^} + {Q} - {f^} - [K] {u^} (18) 

where 

{f^} = (F(u^)} - [K] {u^} + {p^} - Ap(u^) (19) 

[K^] is the nonlinear "reference" matrix, available from the 
elements and Ap(u^) is the change in loads due to grid 
point motions. The advantage to this form is that the second 
order "corrective" load vector {f^} equals zero when the 
forces {F} are linear and when [K] contains linear elastic 
terms (see Eqn 13) . Because the {f} forces are calculated 
at the element level, all known linear elements and super- 
elements may be bypassed in their calculation. This is very 
efficient when a major portion of the structure remains linear. 
However, the advantages for using Eqn (17) are that the load 
error and incremental displacements {u^"**^ - u^} are directly 
available for testing convergence. Therefore, in MSC/NASTRAN, 
both Eqns (16) and (17) are used internally. 

In order to avoid calculating the vector (F(u^)}, 

Eqn (19) is solved for (F(u^)} and substituted into Eqn (12) 
to obtain 

{6^} = {{(u^)} = {p + p’’} - {f^} - [k'^] {u^} (20) 
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Subtracting {6^} for two successive values of i in 
Eqn (20) to eliminate the constants a simple equation is ob- 
tained for the load error in the form 

{5^} = -{f^ - f^'^} - [k] {u^ - u^'^} (21) 

This form provides a simple iteration procedure, namely 

a) For the first iteration after a change in the load vector, 
use Eqn (17) and (20) . Note that the initial values for 
starting from rest are {f°} = {u°} = 0: calculate 

{6°} = {p} + {p*^} and solve [K] {u^} = {6°} (22) (23) 

b) then compute f^ directly from the element and load 
tables : 

= F{u^} - [k""] {u^} + {p'^} - {Ap(u^)} (24) 

c) then compute {6^} from Eqn (21) . 

{6^} = 6° - {f^} + {f°} - [K^] {u^ - u°} (25) 

where {f°} = 0 by assumption. 

d) Use Eqn (17) with i=l again to compute {u } 

[K] {u^ - u^} = 6{u^} (26) 

e) Repeat steps b, c and d to convergence, incrementing the 
superscripts by one at each iteration loop. 

2 . Nonlinear Transient Analysis Using MSC/NASTRAN 

SOL 99 which is a rigid format of MSC/NASTRAN code 
for nonlinear transient problem, is the dynamic comparison to 
nonlinear static systems, providing for combined material and 
geometric nonlinear analysis in the time domain. SOL 66 which 
is a rigid format for nonlinear static problem and SOL 99 
share the same element computer code and solution techniques 
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and provide similar data storage and restart facilities. 

The transient solution is performed in a stepwise manner with 
time steps replacing load steps and with effects of mass, 
damping and unsymmetric matrices added to the stiffness 
matrices . 

Because the calculations are identical, many element 
and material discussions in the previous sections also apply 
to the transient analysis system. The finite elements auto- 
matically provide linear inertial forces and structural 
damping, based upon the initial state, as well as nonlinear 
forces due to geometry and material effects. 

We will develop the basic equations for the transient 
solution. The methods were chosen to be compatible with 
existing MSC/NASTRAN static nonlinear solutions as well as 
with the existing linear transient methods. Thus, the 
Newmark Beta method for transient integration is combined 
with the modifications to Newton's method for nonlinear 
solutions. The additional iteration steps provide equilib- 
rium solutions at each time step, thereby guaranteeing 
stability and accuracy for arbitrary time step size. The 
basic equations are developed below, followed by a discus- 
sion of the stability limits. 

a. Basic Equation 

Adding dynamic loads to the basic equation given 
previously, we obtain a load equilibrium error vector 
at time step n by equation: 
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where 



{6 } = {P - Mu - Bu - F 
n n n n n 



{P„}= 



n n 



{Fn} = 



"Average" load over the time period 

Corresponding acceleration and velocity vectors 
"Average" elasto-plastic element total force 
vector . 



The above equation is solved at the "reduced" 

(u^) displacement vector size. The approximation errors due 
to dynamic reduction methods are not included in the error 
vector Applying Newmark ' s averaging method over a 



finite time period, ^n~l ^ ^n+1 



the "static" forces are 



{F^} = {$F(u^.,) + (l-2B)F(u„) + 3F(u^ ,)} (28) 

n n-rl n n-l 

where 3 = Newmark Beta operator and nonlinear 

force due to a generalized displacement vector . An 

identical definition occurs for from the applied loads 

at each time step. 

From central finite differences, the acceleration 
and velocity vectors are: 



{u } = 
n 



At 



2 = 



At 



o {d } 
2 n 



(29) 



^ = •oTT’ + AU } 

n 2At n+1 n 



(30) 



where {Au ,, ={u ,, -u} 

n+1 n+1 n 

for small angles and At = is the "time step size". 

For large angle changes we require that (u*} , {u}, and {Au} 
are vectors in the global coordinate direction. The vector {u} 
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contains gimbal angles requiring a tranf ormation 
such that 






( 31 ) 



the load error in terms of the current vectors estimates; 

= - [7^2 « " 2^ 'V -[a^ 

At 



{6f„^l + (1-2S)F^ + eF„.^} 



(32) 



and 



{Au , } = {Au + d } 

n+i n n 



(33) 



At any time step the vectors ^^n^ ' ^^n-1^ 

are known. In general the vectors {d }, {u and f ,, = 

n n+i n+i 

F(u^_^j^) must be found either by approximating 
an iterative search. 

Note that for a linear solution {F } equals 

n 

[K]{u^}, where [K] is the reduced stiffness matrix. In this 
case Eqns (32) and (33) may be used to solve for 
directly by setting = 0. With 3=1/3 this will produce 

the linear SOL 27 MSC/NASTRAN results . 

Nonlinear Iterations 

For a nonlinear solution, Equns (32) and (33) may be 
solved with Newton's method (or a modified version) . Using 
{d^} as the primary solution variable provides the following 
iteration algorithm. 
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First assume a linear approximation; 



' {« ) 
n n 



85 



n 



8d 



n 



{d. 



i+1 



n 



- ■= ° 



( 34 ) 



L-t 



{Aj_} = 



85 



n 



8d 



n 



(35) 



Then, the estimated displacement change is 



- d^} = {A, }”-^{5 (d^) } 

n n 1 n n' 



-1 



(36) 



The algorithm is identical to the static case except that 
instead of a tangent stiffness matrix, the "left hand side” 
matrix, obtained from Eqns (32) and (34) is 



{A, }= [-^ M + B + 6K} 



At 



2At 



(37) 



8F 



where [K] = is the current tangent stiffness matrix. 

At a new time step the load iterations may be started by 

O 

assuming that ^^^n+1^ ~ therefore: 



{ 5 ^} = { 0 } 



(38) 



However, to be consistent with displacements, we must extra- 
polate the nonlinear forces; 

= {F„ + (F„ - F„ J } = {2F^ - F^_^} (39) 



n+1 



n 



n n-1 



Substituting Eqns (38) and (39) into Eqn (32) for i=0 results 
in the first estimate; 



{P„} + [AJ {Au„} -{F„} (40) 

n n 2 n n 
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where 



[A^] 




(41) 



Here A is a precalculated matrix. 

previous results. For subsequent iterations, (i>0) Eqn (32) 
becomes : 



{sM = {6°} - 

n n 



[A^Hdi^,} - 



S'^n+1- 






(42) 



where 



[A3] 



M + 
At 




(43) 



At each time step the code will iterate on Eqns 
(36) and (42) until passes the convergence tests or 

the number of passes reaches an iteration limit. With a 
single step, i=l, calculating only 6^, the results will be 
identical to the existing NOLIN results in MSC/NASTRAN. For 
faster convergence the iterations may continue, the matrices 
may be updated, and/or the time step size may be reduced. 

3 . Factors Determining Response to Transient Loading 
The response of a metal to dynamic loading depends 
upon three elements which are characteristic of the material; 
the dynamic constitutive relation, the pressure-temperature 
phase diagram, and the dynamic fracture criteria. In the 
absence of phase changes or fracture, the material response 
is governed completely by the constitutive relation, or 
dynamic equation of state. The constitutive relation is 
characterized as "elastic”, "viscoelastic", "visco-elastic- 
plastic", etc. It depends on macroscopic parameters such as 
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Young's modulus and the viscosity coefficient which, on a 
microscopic scale, arise from atomic interactions and dislo- 
cation processes. The phase diagram determines whether a 
given loading history will bring the metal into a pressure- 
temperature regime where phase transitions are possible. If 
this occurs during shock loading, multiple shock waves may 
result. Microcracks and/or voids begin to appear in a metal 
under tension whenever the criteria for nucleation and growth 
of damage are met. These criteria are functions of material 
parameters and the stress history. 

In this chapter the dynamic constitutive relations 
are discussed. It is the most important element to determine 
the response of a metal to dynamic loadings. 

In macroscopic or rheological language, iron can be 
described as a viscoelastic-plastic material, i.e., it 
responds elastically until yielding (with upper and lower 
yield points) , and is strongly rate-dependent. 

The constitutive relation, or dynamic equation of 
state, thus gives the stress as a function of strain and 

time, or, alternatively, can be written as a differential 

• • 

equation in a,e,a,e and an important goal of present day 
research is to develop the capability to derive this macro- 
scopic constitutive relation by averaging the effects of 
microscopic processes. 

This is analogous to the way in which the thermo- 
dynamic equation of state of a gas is derived by averaging 
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the effects of the molecular collisions. Once developed, 

this capability would be an extremely powerful tool, since 

the response of a metal to dynamic loads could be predicted 

from known microscopic material properties. Although this 

goal has not been attained as yet, there has been significant 

progress during the last few years in this direction for the 

case of crystalline materials, in which dislocation processes 

play the role of the "molecular collisions". In iron, the 

8 2 

density of mobile dislocations is of the order of 10 /cm , 
so the number of dislocations present is certainly large 
enough to make a statistical approach valid. 

The basic approach is based on the relation from 
dislocation theory: [Ref. 3] 

= N bv (44) 

where y is the plastic shear strain rate 
1 ^^i 

(y. . = — t ) where C. is the displacement in 

'll 2 9x . 8x. ' ^ 

D 1 

the x^ direction) , is the mobile disclocation density, 

b is the Burgers vector, and v is the average disclocation 
velocity. This equation provides the link between microscopic 
dislocation motion and macroscopic plastic straint rate, and 
is based on the fact that it has been found possible for many 

p 

crystalline solids to express and v as functions of y 

and applied macrosopic shear stress, x. 
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The constitutive relation is obtained by combining 



Eqn (44) with the incremental Hooke's law 



da_. ^ = 2 ijde?^ + <S 



ID 



ID 



ID 



(45) 



under the assumption that the total strain is the sum of 



elastic and plastic components; 

e P 

e . . = e , . ^ e . . 

ID ij + ID 



(46) 



are the Lame constants, 6 . . is the 

ID 



where A. and p 
Kronecker delta, and a repeated index indicates summation. 

The plastic strain is assumed to cause no volume change, i.e., 
P 



^KK = ® 



(47) 



The combination of Eqns (46) -(48) yields, after taking the 
time derivative : 

a. . -2 u£. • - . - 2 (48) 

13 ^iD KK 13 ID 

In order to combine Eqn (48) with Eqn (44) , we must 
first specify the loading geometry. In plate-slap tests 
under conditions of uniaxial strain, all total strain com- 
ponents except e in the direction of wave propagation are 
zero, and e . Then equation (48) becomes 

Ki\ X 



a - (A+ 2 y) e = -2 y e 

^ X X 



(49) 



Also, in this case, the plane of maximum shear stress in 
isotropic media is inclined at 45*^ to the direction of shock 



propagation, and 



; = i - e P) = 1 

^P 2 ^^x ^ 4 ^x 



( 50 ) 
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The dynamic constitutive relation for this loading geometry 
is thus obtained by combining Eqn (44) , (49) , and (50) : 

If, X, \i , and the right-hand side of Eqn (51) are Icnown 
as functions of macroscopic stress and strain, then the 
equation can, in principle, be solved for any problem 
involving uniaxial strain. 

4 . Dynamic Plastic Deformation of Plates 

We consider cases of the impulsive loading of thin 
metal plates, a simple situation in which an exponential 
shock wave impinges on a clamped square plate. In this case, 
the plastic deformation is supposed to be brought about by 
bending only and combined bending and tension is not analy- 
tically considered. The pressure in the shock wave may be 
assumed to be given by [Ref. 4]. 

P = Pm exp(-t/x) 

where p is the peak pressure at a point at the head of the 

wave, t denotes time and x is a time constant the time 

taken for the pressure to fall to which for high 

explosives is ~30 ysecs. We obtain simply, some of the 
results of Cox and Morland [Ref. 3 ], who theoretically 
investigated the dynamic deformation of plates, side length 
2a, position-fixed at their periphery, when subject to a 
uniform constant transverse pressure p for time x. Their 
detailed examination showed that two categories of load may 
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be distinguished; medium load when p < p < 2p and high 

s s 



load when 2p < p, where p denotes the pressure required 

s s 

just to enforce plastic collapse of the plate, statically 
For p we have, see Fig. 3-1. 



For simplicity we deal only with their first case, 
and for the second, more complicated one, the reader is 
referred to the original source. For this medium range of 
loading two sub-divisions are required: (1) when 0<t<x 

and (2) when t>x. Use is made of the results obtained 
above for static situations using hinge lines, and we take 
over whatever results may be useful. As previously, the 
analysis pays no special attention to material rate effects. 
Further, it is implicit in this approach that the static 
and dynamic modes of deformation are identical. 




s 



(52) 



1 .e . 



Pg = 6Mp/a^ 



Q 

< 

O 




NET ON THE 
FRAME 



TIME 



TIME 



Figure 3.1. Dynamic Loading of Square Plate. 
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Medium Load; 0 < t < x 



Assume the mode of deformation of the plate is 

that of Fig. 3.2 and that the plate diagonals become hinge 

lines. Since p is constant and greater than p , and 

s 

because the four equal parts of the plate act as rigid bodies 

each rotating about its boundary or side, they must undergo 

an angular acceleration, a. If changes in plate geometry 

throughout the deformation are neglected, the resistance to 

deformation, which arises from the yielding taking place in 

the hinge lines, will be constant; this resistance is related 

to p^, the pressure to cause static yield. Thus the excess 

load is constant and hence a must be constant. The work 

done by the excess load will manifest itself as kinetic 

energy acquired by the four rigied regions of the plate, i.e. 
1 ? 

4 * ^ Igj , where I is the moment of inertia of a triangle 
about a side of the plate and o) is its current angular 
velocity. Because a is constant, at the end of time the 
plate centre will have descended through a distance vt/ 2 , 
where v is its current speed. Thus, an equation involving 
the work available for giving rise to rotational kinetic 
energy of the plate triangles is. 



(p - P 3 ) 4a 



2 . ^"^/2 



= 4 • i 



(53) 



Hence substituting for p from (52), noting that I = ma / 6 

s 



and simplifying, 

/ ^^P \ A. niV/.^ ^ (p-Pg) t 

(p - — t = /2 or V = 2 

2 m 



(54) 
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or directly by momentum consideration, thus, 



( (p-p^ ) 



a*2ava _ ma v 
^*2^3 6 ‘ a 



( 55 ) 



The deflection at the end of time t is, 

w, = /2 = (p - p )T^/m 



(56) 



(2) Medium Load: x > t 

After the removal of pressure p at t = x, further deflec- 
tion, w^/ is added to because the rotational kinetic 

energy of the four plate segments must be dissipated in doing 
plastic work in the hinges. Let the segments rotate through 
further angle before coming to rest, then. 



1 w / 

(^ I • — = 4 a/2 • M ((j)>/2) = 8aM • 



Hence 



2„2 

_ ma V 
^2 2 4 M 



mv 



4P< 



/ V 2 2 

<P-Ps> 

m p, 



(57) 



p s s 

using Eqn (55). Thus the total central deflection is. 



2 

w = = p (p - Pg)^ /mPg (58) 

If the total response time of the plate is x^, note that 
x^p^ = pi* In the case of thin plates clamped at their outer 
periphery it is clear that after a significant amount of 
deflection there must be stretching. An analysis of the 
dishing of peripherally clamped metal plates when given a 
transverse initial speed of a few hundred feet per second 
has been given by Hudson [Ref. 3]. Descriptively, Hudson's 
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model is easily understood by reference to Fig. 3-2. Central 
portions of the plate maintain their given transverse speed 
until overtaken by a bending wave propagated radially inwards 
from the clamping ring. Fig. 3-2 (a) shows the imagined plate 
configuration at successive intervals during the deformation 
process; when the plate deflection is significant, the process 
is predominantly one of stretching and for many practical 
purposes the technological approach is adequate. 

Figure 3-2 (a) shows the terminal state taken up 
by a thin rectangular steel plate clamped around its periphery 
after being subjected to the uniform but intense impulsive 
loading, the maximum speed of the plate is a few hundred 
feet per second, but contracting in area with time; the 
dynamic situation during deformation, with plastic hinges 
moving into the plate from the clamped edges and intersecting 
along lines equally inclined to the sides from the corners, 
is indicated diagrammatically in Fig. 3-2 (b) . 

Now the yield criteria of metals should be con- 
sidered. Two approaches to the interpretation of a yield 
criterion are possible; one is Tresca's criterion, the other 
is Von Mises criterion. But for most metals it is generally 
found that while the experimental points fall between the two 
ellipses, they incline towards Von Mises' ellipse. Experiments 
were earlier performed by Lode (1925) on tubes subject to 
tension and internal pressure. Bending and torsion was used 
by Siebel in 1953. All results confirm the general opinion 
that the best simple yield criterion for metals is that of 
Von Mises [Ref. 5]. 
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Figure 3 . 
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2. Dynamic Deformation of Fixed Square Plate. 
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B. PROBLEM DESCRIPTION AND NASTRAN DECK SETUP 



1 . Problem Description 

The model is shown in Fig. 3.3. It is a 20 in. by 
20 in. steel square plate model which is divided to 25 QUAD4 
(NASTRAN card) elements and its boundary is clamped. Its 
thickness is 0.8 inches and the engineering data as follows. 

E = 3.0E+7 
V = 0.3 
L = 20.0 in. 
p = 7.33E-4 Ib/in^ 
t = 0.8 in . 

where E is the elastic modulus, v is the Poisson's ratio, 

L is the length of the plate, p is the mass density, and 
t is a thickness of the plate. 

Next the shock wave applied impinges on the steel 
plate in the -Z direction at time 0.0 seconds (Fig. 3.4) and 
its peak overpressure is 3.0E+4 psi and the shock wave dura- 
tion time (t^) is 5.0E-3 secs. To find the pressure, the 
following formula is used [Ref. 6] 

P ^ exp (-at/ 1^) 

where p^ is peak reflected overpressure at time 0.0 secs 
and a is wave form parameter. But the shock wave pressure 
cannot be applied as a load to the plate since the front face 
shock pressure of the plate is relieved by the rear face shock 
wave which is positive Z direction. The following formula is 
used [Ref . 4 ] . 
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212 , 



p = p exp(-t/T) 

The load data as follows. 

peak lateral loads (-Z direction) 
p = 3 .E+4 psi on the Element (213) 
p = 2.2E+4 psi on the Element (207, 208, 209, 

214, 217, 218, 219) 

Von Mises criterion would be applied for this model. The 
time step (t) to be applied is l.OE-6 secs from 0.0 sec to 
7.0E-4 secs. 
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Figure 3.4 Shock Wave Applied. 
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Figure 3.5. 
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2. NASTRAN Input Data Deck 



Reference is made to the echo of the input data deck 
which appears in the output. 

a. Sol 66 [Ref. 6] 

1 NASTRAN PREF0PT=2 

2 ID NONLINEAR, STATIC 

3 SOL 66 

4 TIME 30 

5 CEND 

6 TITLE-NONLINEAR DEFORMATION WITH STATIC LOADS 

7 SUBTITLE=TRANSIENT RESPONSE BY SHOCK WAVE LOAD. 

8 ECHO=UNSORT 

10 DISP=ALL 

11 OLOAD=ALL 

12 SPCF=ALL 

13 GPFOR=ALL 

14 ELST=ALL 

15 SEALL=ALL 

16 SPC=200 

17 SUBCASE 1 

18 LOAD=10 

19 NLPARM=10 

20 BEGIN BULK 

21 $ 

22 $ LATERALLY LOADED FLAT PLATE 

23 $ 
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24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 
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MESHOPT, , , , ,YES 

$ 

$ USE MSGMESH TO GENERATE 5*5 ARRAY OF QUAD4 ELEMENTS 

$ 

EGRID, 1, , 0 . , 0 . 

EGRID,2, ,20. ,0 . 

EGRID,3, ,20. ,20. 

EGRID,4 , ,0 . ,20. 

GRIDG , 2 , , , 5 , -1 , -2 , -3 , , +FLD2 
+FLD2, 5 , -4 

CGEN,QUAD4 , 201,200 , 2 

$ REFER TO PAGE 2 . 18-7 , APPLICATION MANUAL. 

PSHELL,200 ,65 ,0 . 8, 65, ,65 
MATl,65,3.0E+7, , .3,7.33E-4 
MATS1,65,40,PLASTIC, ,l,l,2.6E+4 
TABLES 1,40 ,0 ,,,,,, ,+TABl 

+TAB1 ,0.0,0.0,0.001,3. OE+4 ,0.025,4. -E+4 ,0.05,5. 8E+4 , +TAB2 

+TAB2 ,0.10,6. 7E+4 ,0.125,7. OE+4 ,0.2,7. 5E+4 , ENDT 

PARAM , COUP MASS , 1 

PARAM,K6ROT,1.0 

PARAM,LGDISP,1 

PLOAD2, 10, -14. 7, 208, 212, 213, 214, 218 
NLPARM,10 ,1, ,AUTO, , 20 ,W,YES 
$NLPARM,20,3, ,AUTO, ,20,W,YES 
SPCG,200 ,2,123456,A,B 
SPCG,200 ,2, 123456,B,C 
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SPCG,200,2,123456,C,D 

51 SPCG,200,2,123456,A,D 

52 ENDDATA 

1) Executive control card 

card 1 - The NASTRAN card contains the keyword PREFOPT 
= 2 which requests that MSGMESH be used to 
generate input data. 

card 3 - Request nonlinear static analysis, Rigid format 

66 . 

2) Case control card 

card 15 - Requests superelement generation and assembly 

* should be set as 'all' 

card 18 - Requests static load combination. 

* SID should be matched to the PLOAD2 card in 
the bulk data deck. 

card 19 - Requests nonlinear parameter selection. 

* SID should be matched to the NLPARM card 
in the bulk data deck. 

* This card may appear above or within a 
subcase. 

* LOAD card is only applicable in statics, 
buckling, heat transfer problems. 

3) Bulk data cards 

card 24 - Requests out options control 

* An entry of YES causes the printing of the 
' input card echo ' . 
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card 28-31 - Defines vertices of quadrilateral field 

using MSGMESH EGRID command, 
card 33 - Defines a field of grid points and generate 
grid cards. 

* The grid number generated by MSGMESH starts 
from the FID of GRIDG (if FID is 2, the first 
grid point ID is 20000) 

* 5 means 5 elements in each direction, 
card 34 - Generate element connection. 

* QUAD4 is the element type which should be 
generated. 

* Element ID number starts from 201. 

* Fifth field 200 is the PID of pshell card, 
card 38 - Specifies table references and material 

properties which are stress-dependent for use 
in material nonlinearity applications 

* Taken Von Mises yield function, isotropic 
hardening rule. 

card 39 - Represent stress-strain data. 

card 43 - The value 1.0 suppresses singularities 

* It is intended primarily for geometric 
nonlinear analysis. 

card 44 - The value 1 allows the nonlinear calculations, 
card 45 - Request the initial condition on the elements 
in the negative Z direction, 
card 46 - Controls the nonlinear analysis iteration. 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 



* To reduce the output, use sixth field. 

* For rapid convergence, increase the Error 
tolerance and maximum iteration and the limit 



on diverging iterations, 
card 48-51 - Generate spcl cards. 

This card is only used to MSGMESH generation, 
b. Sol 99 
NASTRAN PREFOPT=2 
ID NONLINEAR, TRANSIENT 
SOL 99 
TIME 59 
CEND 

TITLE=N0NLINEAR DEFORMATION WITH DYNAMIC LOADS 
SUBTITLE=TRANSIENT RESPONSE BY SHOCK WAVE LOAD. 
ECHO=SORT 

SET 1=20202,20203,20204,20304,20404 

SET 2=201,208,212,213,214,219,225 

SET 3=20000,20003,20005,20200,20400 

SDISP=1 

OLOAD=l 

SPCF=3 

ELST=2 

SEALL=ALL 

SPC=200 

LOADSET=77 

SUBCASE 1 
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20 

21 

22 

23 

24 

25 

26 . 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 



DLOAD=37 



TSTEPNL=45 
SUBCASE 2 
DLOAD=37 
TSTEPNL=46 
SUBCASE 3 
DLOAD=37 
TSTEPNL=47 
OUTPUT (XYPLOT) 

PLOTTER NAST 
CSCALE 1,3 
XTITLE=TIME IN SECS 
XGRID LINES=YES 
YGRID LINES=YES 

YTITLE= DISPLACEMENT GRID 20202 
XYPLOT SDISP RESP/20202 (T3) 

YTITLE=VELOCITY GRID 20202 
XYPLOT SVELO RESP/20202 (T3) 

YTITLE=DISPLACEMENT GRID 20101 20102 20202 
XYPLOT SDISP RESP/2010KT3) ,20102(T3) ,20202 (T3) 
BEGIN BULK 
$ 

$ LATERALLY LOADED FLAT PLATE 

$ 

MESHOPT, , , , ,YES 

$ 
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46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 



$ USE MSGMESH TO GENERATE 5*5 ARRAY OF QUAD4 ELEMENTS 

$ 

EGRXD /0» /0» 

EGRXD /2^ /20» /0» 

EGRXD, 3, ,20. ,20 . 

EGRXD, 4, ,0. ,20. 

GRXDG ,2, , ,5,~1,~2,~3, , +FLD2 
+FLD2,5,-4 

CGEN,QUAD4,201,200,2 

$ REFER TO PAGE 2 . 18-7 , APPLXCATXON MANUAL 
PSHELL,200,65,0.8,65, ,65 
MAT1,65, 3.0E+7, , .3,7.33E-4, , ,0.01 
MATS1,65,40,PLASTXC, ,l,l,2.6E+4 
TABLESl , 40,0,,,,,,, +TAB1 

+TAB1 ,0.0,0.0,0.001,3. OE+4 ,0.025,4. OE+4 ,0.05,5. 8E+4 , +TAB2 

+TA32 ,0.10,6. 7E+4 ,0.125,7. OE+4 ,0.2,7. 5E+4 , ENDT 

PARAM , COUPMASS , 1 

$PARAM,AUTOSPC,YES 

PARAM, W4,1.6l6E+3 

PARAM, K6 ROT, 1.0 

PARAM, LGDXSP,! 

PARAM, SLOOPXD,! 

P ARAM, LOOP XD,1 
PARAM, STXME,0 .000450 
PARAM,MAXLP,7 

DLOAD ,37,1.0,1.0,7,1.0,6,1.0,5 
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72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 



TL0AD1,7,61, ,0,22 
LSEQ,77,61, 71 
PLOAD2 ,71,3. 2E+4 , 213 
TL0AD1,6,62, ,0,22 
LSEQ, 77,62, 81 

PLOAD,81,2.5E+4,207,208,209,212,214 
TL0AD1,5, 63, ,0,22 
LSEQ, 77, 63, 91 

PLOAD2,91,2. 5E+4,217,218,219 
TABLEDl , 22 ,,,,,,, , +NEXT1 

+NEXT1, -1. ,0.0,0.0,-0.0010,0.0001,-1.0,0.0010,-0.560, +NEXT2 
+NEXT2,0 .004 ,-0 . 001,ENDT 

TSTEPNL,45,150 ,1.0E-6,20,AUTO, ,20 ,W,+LLl 

+LL1 , 1 . OE-7 , 1 . OE-7 , 1 . OE-7 ,20,20,10 

TSTEPNL ,46,200,1.0E-6,20, AUTO , , 20 , W , +LL2 

+LL2,1.0E-5,1.0E-5,1.0E-5,20,20, 10 

TSTEPNL , 4 7 , 120 , 3 . OE-6 , 20 , AUTO , , 20 , W , +LL2 

+LL2 , 1. OE-4 , 1. OE-4 , 1. OE-4 ,20,20,10 

SPCG,200,2,123456,A,B 

SPCG, 200,2,123456,B,C 

SPCG,200,2,123456,C,D 

SPCG, 200 ,2,123456 , A, D 

ENDDATA 

1) Executive control cards 

card 3 : Requests nonlinear transient analysis. Rigid 

format 99. 
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2) Case control cards 

card 18: Selects a static load set for use in dynamics. 

Thus this may be referenced by dynamic loads, 
card 20 : This card selects a dynamic load for nonlinear 

transient problem. 

card 22-24: Requests restart at time 0.00015 secs, 
card 25-27: Requests restart at time 0.00035 secs, 
card 28-39: Requests the xy plotting which is 

corresponding to the time step. 

3) Bulk data cards 

card 64 : This gives damping effects for transient 

analysis . 

* Sometimes, more rapid convergence effects 
can be obtained by this damping factor with 
damping coefficient in MATl card (GE) . 

card 67: This card identifies the initial conditions 

from a previous SOL 66 nonlinear static 
solution. 

* Refer to the param card, SLOOPID. 

card 68,69: Request these two cards when we taken more 

than one subcase . 

card 70: Requests more internal calculation. 

* We can reduce iteration steps and prevents 
the overprinting due to iteration steps. 
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card 71: 



Defines a dynamic loading condition for 
transient response problems as a linear 
combination of load sets 



* See the remarks 5, 7, 9 of DLOAD card in 
bulk data. 

card 72, 75, 78: Defines time dependent dynamic load. 

* Requests a static load card set (LESQ) . 

card 73, 76, 79: Defines a sequence of static load sets. 

This load sets may be referenced by dynamic 
load card (TLOADl) 

* This card should be used when we want to 
apply area pressure on elements than a point 
pressure on grid point. 

card 74, 77, 80: Defines a uniform static pressure load 

on the plate elements. 

card 81-83: Defines a tabular function for use in 

generating a time-dependent dynamic loads. 

* In this card the peak pressure (3.0E+4 psi) 
occurred at O.OOOl sec. 

card 84, 86, 88: This card control the nonlinear 

transient analysis. 

* This card is expalined in the sensitivity 
analysis . 

C. SENSITIVITY ANALYSIS 
1. General 

In transient analysis, two types of instability could 

occur. The first is the familiar nonlinear load iteration 

divergence which also occurs in static analysis. The second 
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is the divergence which grows with time in the transient 
integration. Both instabilities are caused by uncorrected 
nonlinear equilibrium errors . 

A convenient method, from Von Neumann, for analyzing 
the stability limits is to assume that the nonlinear forces 
are nearly linear and the error vector has a constant conver- 
gence. It is assumed that the nonlinearity has a first order 
approximation; [Ref. 7] 

{F(y)} '' [K^ + {y} (59) 

NL 

where AK is the difference between the tangent stiffness 
matrix and its approximation K . The error vectors, {6}, 
are assumed to grow at the rate A, defined as: 

{5^'^^}= X{6^} (60) 

Note that if lX| >1 the system will be defined as unstable. 

After lengthy calculations we may summarize the various 
criteria, assuming that the matrices are reduced to equivalent 
scalar modal quantities, the criterian for stable solutions 
are ; 

1) For time steps with converged "static" iterations: 

1/2 > 3 > 1/4 

This is the same criteria as linear analysis. 

2) For "static load" iterations; 

3 A < 3 B + M 

At^ 

This states that the mass and damping add to the effective 
linear stiffness [K^] and improve the stability. 
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3) For time step integration with no intermediate static 
iterations, as in the standard Newmark Beta method: 

NL r 4 

A K < (4 8 -1) M 

At 

This restriction is more severe than criteria b above, 
proving that the internal iterations are more stable 
than the Newmark integration. 

In summary, the Von Neumann method will have fewer 
divergence problems than either the static nonlinear solu- 
tion or the single step transient nonlinear methods. The 
better stability of the method and the capability to use 
larger time steps outweighs the cost of a few internal 
iterations on the static element forces. 

2 . Nonlinear Transient Controls 

The input data required for SOL 99 is a combination 
of transient control data, similar to SOL 69 (Direct Transient, 
superelements) , and nonlinear modeling data similar to SOL 66 
(Nonlinear Statics) . The nonlinear properties are defined by 
nonlinear material data (MATSl) , gap elements (GAP) , and 
large displacement effects (parameter LGDISP) . The transient 
effects are produced by loading functions (TLOADi, DAREA, etc.) , 
damping (parameters, elements, and material data) , and element 
mass properties. 

The only unique data required for SOL 99 is supplied 
on the TSTEPNL data cards. (Several optional parameters are 
provided for tuning the method and for use in restarts.) The 
TSTEPNL card in itself is an amalgam of the normal transient 

TSTEP card and the nonlinear control card, NLPARM. 
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a. Case Control 



The significant rules for nonlinear transient are: 

1) Each subcase defines a time interval starting from the 
last time step of the previous subcase which is subdivided 
into small time steps. The output TIME printout will be 
labeled by the net time, including all previous subcases. 

2) The input loading functions may be changed for each sub- 
case or continued by repeating the same DLOAD request. 

The bulk functions are evaluated at the net model time, 
which starts at the last time step of the previous subcase. 

3) Each subcase may have a different time step size, time 

interval, and iteration control selected by the TSTEPNL 
request. The case control requests which may not be 
changed after the first subcase are: SPC, MPC , DMIG, and 

TF. 

4) Output requests for each subcase are processed indepen- 
dently and the output for each case are generated before 
the next subcase is processed. Separate XYPLOT ouputs will 
occur for each subcase. Available outputs are DISPLACEMENT, 
VELOCITY, ACCELERATION, OLOAD , STRESS, FORCE, S DIS PLACEMENT , 
SVELOCITY, and, SACCELERATION . Note that NONLINEAR and 
SPCFORCE outputs are not provided. 

5) Initial conditions are not explicitly provided for user 
input in SOL 99. The reason is that the nonlinear element 
configuration and equilibrium state must be known at every 
time step. If initial conditions were given, all of the 
nonlinear element forces and plastic stresses, as well as 
the past history, must be given. Therefore, the initial 
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conditions become part of the analysis and are generated 
by ini.tial transient subcase or by restarts from a previous 
static solution. 

6) All normal superelement model generation options and 

matrix reduction options are allowed for the linear portion 
of the structure. General Dynamic Reduction, Component 
Mode Synthesis, and Guyan Reduction may be performed for 
nonresidual superelements. Therefore, the residual super- 
element may contain scalar degrees of freedom representing 
linear model formulations, 
b. TSTEPNL Data 

The TSTEPNL data card description is provided in 
Section 2 of the MSC/NASTRAN User's Manual. The input fields 
specify the time step size, the number of steps, and the out- 
put interval as well as the nonlinear iteration options. 

The nonlinear iteration data are provided for 
controlling the accuracy and stability of the nonlinear 
calculations and they will not effect linear solutions. The 
controls are similar to those on the NLPARM bulk data used in 
static nonlinear analysis. 

As discussed above two types of errors or insta- 
bility are encountered in the transient solution. The load 
iterations are similar to static load analysis and may diverge 
at one or more time steps when large nonlinear stiffness 
changes occur. Load iteration divergence may be cured by 
smaller time steps or by stiffness matrix updates. The 
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diagnostic output value LAMBDA-I (from DIAG 50) indicates 
the rate of growth of this error. 

The second type of divergence occurs in the time 
domain when the load iterations do not converge, or are limi- 
ted by the user. The errors may increase for each new time 
step until they cause the code to abort. This problem may be 
cured by allowing more load iterations or by using smaller 
time steps. The diagnostic output value LAMBDA-T (from DIAG 
50) indicates the rate of growth of this error. 

The size of the time step is of primary importance 
in transient analysis to control -accuracy and avoid diverging 
solutions. The use of several subcases with different time 
step sizes to correspond with nonlinear activity is recommended 
as a normal procedure. When the user wishes to avoid overly- 
small time steps for highly nonlinear zones, he may increase 
the iteration limits and/or the error tolerance and the diver- 
gence limits. 

c. Restarts and Data Base 

Solutions 66 and 99 share the same data base/ 
superelement file storage logic for the nonlinear tables and 
matrices. Therefore, the restart system for transient analysis 
is allowed to use either a previous static or transient non- 
linear analysis as initial conditions. 

For restarting from previous static analysis only 
the first subcase is affected. Simply provide a data base 
created in SOL 66 and specify the parameter. 



65 



PARAM,SLOPPID,N 



where N was the printed value of LOOP for the desired 
static solution. Constraint sets, direct input matrices, 
mass, and damping may be changed. 

Restarts from a previous transient execution are 
available for a number of cases. If the same model is to be 
reexecuted, only the residual superelement needs to be re- 
assembled (SEMA, SELA=0) . If the end results for the previous 
transient run are to be used for the initial conditions at 
t=0 , add the parameter LOOPID=N where N corresponds to the 
Nth subcase of the previous run and N dummy SUBCASE commands 
to start the residual case control execution. 

The normal restart for a transient run is to 
continue from the last step of a previous subcase with 
diff erent . loads and/or TSTEPNL data. Provide the following 
parameters : 

LOOPID = N - Start from the Nth subcase 

Stime = t - Start from time t. 

Note that constraint sets should not be changed 
to avoid incompatible matrix sizes. 

D. NUMERICAL RESULTS 

The progressive deformations of the plate are plotted 
along the x-axis at five different time locations as shown 
in Fig. 3.6(b). At the early time (1.5E-4 secs), the response 
of the plate is in the elastic region. During the deformation, 
plastic hinge lines are formed along the clamped edges and the 
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diagonals of the plate as indicated in Fig. 3.6(a). This 
response phenomenon agrees well with Hudson's studies 
[Ref. 4]. As described in Section A-4 of Chapter 3, when 
the plate deflection is significant, the plate experiences 
in-plane stretching. 

Time history displacement responses of Grid points 
20101, 20201, and 20202 are plotted in Fig. 3.7. As shown 
in the figure, for the short duration impulse, the maximum 
displacement occurs at the center of the plate after termi- 
nation of shock loading. The rate of increase in displacement 
is the maximum at the termination of shock loading' as shown 
in Fig. 3.8. During the application of the shock loading, 
the acceleration level is high, implying the high inertia 
force in the system. At the termination of the inertia 
force in the system, the acceleration stays in the constant 
level until the plate displacement reaches the maximum. 

Time hisotry Von-Mises stress of QUAD4 element 213 (plate 
center) is plotted in Fig. 3.7. The stress variation is 
highly non-linear due to the material non-linearity. 
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displacement (Inches ) 




Figure 3.6(a). Plastic Hinge Line Formation. 
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Figure 3.6(b). Progressive Deformation of the Plate 

Along Section A-A. 
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Figure 3.7. Time History Displacement Responses of Grid Points 

20101 , 20201 , 20202 . 
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Figure 3.8. Time History Velocity Responses of Grid Point 20202. 
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Figure 3.9. Time History Von-Mises Stress Response of QUAD4 Element 213 



IV. CONCLUSIONS AND RECOMMENDATIONS 



The state-of-the-art literature review revealed that the 
skirt plate is very effective to absorb initial explosion 
energy before the shock loading reaches the main armor. The 
material of the armor must possess the high hardness to prevent 
the severe damage and the high ductility to minimize the 
spalling. The composite materials may serve more effectively 
than the conventional material to minimize the damages against 
the powerful modern projectiles. 

The non-linear finite element analysis of the plate sub- 
jected to the air shock loading verifies the plastic hinge 
line formation and in-plane stretching inducing snap-through 
collapse. For the short duration shock loading, the maximum 
deformation occurs after the termination of loading, and the 
velocity reaches the maximum at the termination of shock. 

The Von-Mises stress is highly non-linear due to the material 
non-linearity. For the non-linear analysis using MSC/NASTRAN, 
TSTEPNL bulk data card is found to be the most important card 
which controls the accuracy and the stability of non-linear 
calculation . 

To insure the convergence of solutions, convergence limit 
in the TSTEPNL card in the bulk data deck should be considered. 
If a program was run with big convergence limit, the solution 
will be gotten rapidly with small iteration steps but it is 
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not accurate and large time steps can't be taken. On the 
contrary, if the program was run with small convergence 
limit, the solution will be accurate but it is gotten slowly 
with large iteration steps. Thus, it is expensive and anyway 
large time steps can't be taken too. Thus, compromise should 
be taken to run a large time steps. For example, if one non- 
linear transient problem is needed total 700 time steps to get 
the solutions, the following procedure is recommended. In the 
first subcase, the program is iterated with smaller conver- 
gence limit, lO*E-7 to get more accurate solution and in the 
second subcase, the program is run with limit, 10*E-5, and in 
the third subcase, the program is iterated with convergence 
limit, 10*E-4. 

Consequently, the total 700 time step will be calculated 
successively . 

The sensitivity studies on shock period and amplitude is 
suggested to understand their effects on plate responses. 
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APPENDIX A 



THE COMMAND TO OPERATE DATA BASE 



****THE COMMAND TO ERASE THE DATA BASE NAMED "MSS . S2 770P . NODE " 

//ERASE JOB (2770,0252) ,NASTRAN,CLASS=A 
//*MAIN 0RG=NPGVM1.2770P 
// EXEC PGM=IEFER14 
//SYSPRINT DD SYSOUT=A 

//DDl DD DSN=MSS .S2770 .NODE, DISP=(OLD, DELETE) 

// 

****THE COMMAND TO PLOT FROM THE DATABASE NAMED "MSS . S2 770 . PPLOT" 

//LEEPLCT JOB ( 2 770 , 02 52 ) , NASTRAN , CLASS=C 

//*MAIN ORG=NPGVMl . 2 77 OP 

//♦FORMAT PR, DDNAME=SYSVECTR, DEST=LOCAL 
// EXEC NASTPLOT,PLTDSN= 'MSS .S2770 .PPLOT' 



74 



appendix b 



DATA DECK FOR PROGRAM 



general description of data oecx 

The Input deck begins with the required resident operating system control cards. The type 
and ninber of these cards will var 7 with the Installation. Instructions for the preparation of 
these control cards are given In Section 7.6 of the MSC/NASTRAH Application Manual. 

The operating system control cards are followed by the MSC/NASTRAN Data Deck, which consists 
of the following three sections: 

1. Ejtecutive Control Deck 

2. Case Control Deck 

3. Bulk Data Deck 

The E;«ecut1ve Control Deck and the Case Control Deck both have free-fleld formats- Only 
coluins 1 through 72 are used for data. Any Information In columns 73 through 30 may appear In 
the printed echo, but the data will not be used by the program. As explained In Section 2.A.1, 
limited use Is made of the data In columns 73 through 30 for the Bulk Data Deck. If the last 
character on a card Is a comma (not necessarily In column 72), the next card Is a continuation of 
this physical card. Any number of continuation cards may be specified, and together they form a 
logical card . 

The NASTRAH card Is used to change the default values for certain operational parameters, 
such as buffer size or the number of data lines printed per page. More than one NASTRAN card may 
be present. NASTRAN cards are optional, but. If present, they must be the first cards of the data 
deck. The NASTRAN card Is a free-fleld card (similar to cards In the Executive Control Deck). 
Its format Is as follows: 



NASTRAN keywordj^ » value, keyword 2 • value, ... 
An alternate format Is: 



NASTRAN keyword^ « value 
NASTRAN keyword 2 ■ value 
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Input Data Card EG3I0 Grid Point Position 



Description : Defines the position of a geometric grid point of the structural model 



Format and Example: 



1 


2 


J 


4 


5 


6 


7 


a 


9 


10 


EGRID 


ID 


CP 


XI 


X2 


X3 






jXi 


i 


EGRID 


2 


i 3 i 

1 


j 1.0 


2.0 


3.0 











Field 


Contents 






ID 


Grid point Identification number {Integer > 0) 






CP 


Identification number of coordinate system In which the 
Is defined (Integer ^0 or blank) 


position of 


the grid point 


X1.X2.X3 


Position of the grid point In coordinate system CP (Real 


or Blank) 




Remarks: 1. 


The meaning of XI, X2, and X3 depends on the type of 
follows: (see C3RD1 card descriptions) 


coordinate 


system, CP, as 



Type 


XI 


X2 


X3 


Rectangular 


X 


Y 


Z 


Cyl indrical 


R 


9(degrees) 


Z 


Spherical 


R 


9 (degrees) 


t(degrees) 



2. The EGRID card Is used In MSGJ^ESH to define curved edges and for other special 
purposes (see Part I, Section 3.0). It Is similar to a GRID card, except that It 
win not be processed by MSC/MASTRAM. 

3. Entries In fields 7, 8 and 9 are not checked by MSGMESH. 
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Input Data Card GRIOG Linear and Higher Order Interpolation 



Description : Defines a field of grid points and generates GRID cards 



Format and Example : 

1 23456789 10 



GRIDG 


FID 


CD/CI 


PS 


L 


GA 


GB 


GC 


SEID 




GRIOG 


10 


*0 


18 


-1 


-2 


-3 




*G1 




M 


GO 


H 


GE 


GF 


GG 


GH 


Cl 




+G1 


9 


-4 














4-02 




E1(A8) 


£2(AB) 


a(BC) 


£2(BC) 


El(CD) 


E2(CD) 


EKOA) 


E2(DA) 


i 


-t-GZ 










6 


5 


8 


i 7 






EKEF) 


£2(EF) 


Q(FG) 


E2(FG) 


EKGH) 


E2(GH) 


El(HE) 


E2(HE) 
























EKAE) 


E2(AE) 


EK3F) 


E2(BF) 


EKCG) 


£21 CG) 


El ( DH ) 


E2(DH) 

























Field 

FID 

CO/CI 



Contents 

ID number of a GRIDG field (Integer, 1 ^ FID ^ 99). Must be unique for all GRIDG 
cards. When entry is negative, internal EGRID cards are generated. 

Prescribes displacement and/or interpolation coordinate systems (Integer or 
blantc). See Remarks 4 and 5. 



PS 

L.M.N 

SEID 

GA , . . . , GH 



PS entry on generated GRID cards (any of the digits 1-6 with no imbedded blanks) 
( Integer ^ 0 or blank) 

Grid field dimensions when an entry, L, M, and/or H is positive. (See Figure 3.) 
Absolute value is a LID entry on a LIST input card, when an entry is negative. 

ScID entry on generated GRID cards (Integer >_ 0 or blank) 

Grid point identification numbers, GIDs, of grid field vertex points (Integer / 0 or 
blank). See Remarks 3 and 6. 



Cl CP entry on generated GRID cards* (Integer ^ 0 or blank), and grid field inter- 

polation coordinate system definition. See Remarks 4 and 5. 

El(XY) ID number of edge grid point (closer to X) on side XY (Integer A 0 or blank). See 

Remarks 3 and 6. 

E2(XY) ID number of edge grid point (closer to Y) on side XY (Integer A 0 or blank). See 

Remarks 3 and 6. 



(Continued) 
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Remartcs: 1. 

2 . 



See Part I, Sections 2 and 3 for discussions on the GRIDG Input card. 
Required data by type of field 



Type 


Required Data 


LI. 'IE 


FID, L, GA, GB 


TRIA 


■ FID, L, GA, G3, GC 


QUAD 


FID, L, M, GA, GB, GC, GD 


HEX 


FID, L, M, N, GA thru GH 



Only the parent and first continuation card are applicable for linear 
Interpolation. 

3. Vertex and edge points are defined on GRID or EGRID input cards, or on GRID cards 
generated by previously processed GRIDG input cards. Field vertex points are 
automatically equivalenced to corresponding entries GA, ..., GH, when these entries 
are positive. GRID cards are generated only for active field grid points. See 
Part I, Section 10. 1. 

4. CD/CI entry may be used to prescribe both the displacement coordinate system, CD, 
and the Interpolation coordinate system. Cl. Both CD and Cl are prescribed when 
the entry is positively signed. Only CD Is prescribed when the entry is 
unsigned. Only Cl Is prescribed when the entry is negatively signed. The Cl entry 
on the first continuation card may be used when both CD and Cl exist and are 
di fferent. 

5. Cl references a CID entry on a coordinate system definition Input card. When Cl is 

entered as 0, Interpolation is in basic rectangular coordinates. When Cl is not 
entered. Interpolation Is in the coordinate system that positions the field vertex 
and edge points. In this case, all vertex and edge points trust be defined in the 
same coordinate system (C? entries on GRID and/or EGRID cards). See Part I, 
Section 3.1 and Table 2. 

6. The > Interpolation point components, In spherical coordinate systems, or 

the 9 Interpolation point components, in cylindrical coordinate systems, may be 

offset by a i360® by appending a tS or a *C to grid point entries. 

7. If an entry for GD is present In a TRIA field, grid point D replaces grid point C 

In Interpolating grid point positions from side CA. Grid point D is usually 

defined by an EGRID card. M must be blank. See Part I, Section 3.5. 

3. If both El(XY) and E2(XY) are entered, the positions of generated grid points are 

determined by cubic Interpolation from edge XY. If £2(XY) Is blank, positions are 
determined by quadratic Interpolation, If both El(XY) and E2(XY) are blank, 
positions are determined by linear Interpolation. Either none or one edge point, 
E(XY) may be specified for each edge of a TRIA field. El(CD) positions a grid 

point on edge CA of TRIA fields. See Part I, Section 3.7. 

9. Intermediate continuation cards cannot be omitted even if all entries are blank. - 
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Input Data Card CGEM Generate Element Connection Bulk Data Cards, Form 1 
Description ; Generate element connection cards for element types listed under Remark 1 



Format and Examole: 



1 


2 


3 


4 


5 


S 


7 


8 


9 


10 


CGEH 


npE 


FEID 


PID 


FID 


DIR 


TH 


EIDL 


EIDH 




CGEN 


QUA04 


100 


10 


1 










+ C2 






TA 


T3 


TC 


FD 






+C2 


0.002 


0.002 


0.001 


0.001 








1 


1 



Field Contents 

TYPE Any of the mnemonics listed under Remark 1 



FEID 

PIO 



The absolute value of FEID is the element Identification number, EIO, of first 
element In element set (Integer ^ 0). See Remark 4. 

ID of an element property card or PGEN card (Integer > 0 or blank). Must be unique 
with respect to all other property card IDs and PGEM IDs. See Remark 5. 



FID ID of associated grid point field defined on a GRIDG or GRIDU card 

DIR One of the entries L, M, or M (Default entry Is L). Direction of elements In the 

field. See Table 4 of Part I, Section 7.1. 



TH Material prooerty orientation specification (Real or Blank; or Integer ^ 0 and less 

than 1,000,000). If Real or Blank, specifies the material property orientation 

angle In degrees. If Integer, the orientation of the material x-axis is along the 

projection onto the plane of the element of the x-axis of the co-ordinate system 
specified by the Integer value. 



EIDL 



ID of element, below which no element connection cards are generated (Integer > 0 
or blank) 



EIDH 



ID of element, above which no element connection cards are generated (Integer > 0 
or blank) 



TA,TB,TC,TD 



Thickness of corner points of QUAD and TRIA fields when TYPE Is TRIA3, TRIA6, CIIAD4 
OP (JUA08 (Real or Blank). See Remark 6. 



Remarks: 1. 



Element types: R0D, PL0TEL, CFTU8E, TUBE, TRIA3, TRAPRG, TRIARG, 

TRIAX6, (JUAD4, SHEAR, QUADS, HEXAl, HEXA2, HEXA8, HEXA20 and HEX20. 
HEXA8 and HEX^O refer to HEXA elements with 8 or 9 through 20 nodes. 



TRIAS, TRIMS, 
The mnemonics 
respectively . 



(Continued) 
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2. Elements are tabulated by class In Table 3, of Part I, Section 7.0. Line class 
elements may be generated In all field types. Tr1a3 and Tr1a6 class elements may 
be generated In TRIA and (XIAO fields. (Juad4 and Quad3 class elements may be 
generated In TRIA, QUAD and HEX fields. Only QUAD4 and QUAC8 element types are 
valid In TRIA fields. TRIA3 and TRIA6 elements are generated. In this case, along 
the field edge corresponding to the element direction. Hex8 and Hex20 class 
elements may only be generated In HEX fields. 

3. MSGMESH automatically generates PENTA elements In place of HEXA and HEX20 elements, 

when adjacent element edges would have the same grid point Identifiers. TRIA6 

elements replace QUADS elements and TRIA3 elements replace 0UAD4 elements, when 
adjacent element vertices would have the same grid point Identifiers. Otherwise, 
MSG.MESH suooresses the generation of an element which would reference Identical 
grid point Identifiers. Elements are generated only when all vertex points 
exist. Deleted midside nodes suopress the generation fo TRIMS and TRIAX6 elements. 

4. EIOs are Incremented when the FEID entry Is positive and decremented when the FEID 
entry is negative. 

5. The PID entry may be blanic for those elements whose element connection card permit 
a blank property ID. 

6. Element thicknesses may be generated for TRIA3, TRIAS, QUA04 and QUADS elements 
located In TRIA and QUAD fields. Interpolation from the thicknesses at the field 
vertices Is based on field topologies with equal spaclngs of thicknesses for both 
GRIOG and GRIDU fields. Unequal spaclngs of thicknesses are defined for GRIDG 
fields specified with unequal spaclngs of grid points. 
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Input Data Card MATl Material Property Oefinitlon, Form 1 



Description : Defines the material properties for linear, temperature-independent, isotropic 

materia is 



Format and Example; 



1 


2 


3 


4 


5 


6 


7 


a 


9 


10 


MATl 


MID 


C 


G 


W 


RHO 


A 


TREF 


GE 




MATl 


17 


3. *7 




0.33 


4.28 


6.5-6 


5.37^2 


0.23 


ABC 






ST 


SC 


SS 


MCSID 












*0C 


20. f4 


IS. •►A 


12. +4 


1003 













Field 

MID 

E 

G 



NU 



RH9 

A 

TREF 

GE 

ST,SC,SS 



MCSID 



Contents 

Material Identification number (Integer > 0) 

Young’s modulus (Real or blank) 

Shear modulus (Real or blank) 

Poisson's ratio (-1.0 < Real _< 0.5 or blank) 

Mass density (Real) 

Thermal expansion coefficient (Real) 

Thermal expansion reference temperature (Real) 

Structural element damping coefficient (Real) 

Stress limits for tension, compression, and shear (Real). (Used only to compute 
margins of safety in certain elements; they have no effect on the computational 
procedures. ) 

Material Coordinate System identification number (Integer > 0 or blanx) 



Remarks: 1. 

2 . 

3. 

A. 

5. 

6 . 



The material identification number must be unique for all MATl, MATE, MAT3 and MAT9 
cards. 

MATl materials may be made temperature dependent by use of the MATTl card. 

The mass density, RH(J, will be used to automatically compute mass for all 
structural elements. 

'height density may be used in field 5 if the value l/g is entered on the PARAM card 
MTMASS, where g Is the acceleration of gravity (see Section 3.1.5). 

MCSID must be nonzero If the CURV module Is used to calculate stresses or strains 
at grid points on plate and shell elements only. 

To obtain the damping coefficient, GE, irultiply the critical damping ratio C/C„, by 

2 . 0 . 



(Continued) 



81 



MATl (Cont.) 



7. Either E or G fnust be specified (I.e., nonblank). 

8. If any one of E, G, or NU is blank, it «il 1 be computed to satisfy the identity 
£ « 2(1*NU)G; otherwise, values supplied by the user will be used. This 
calculation is only oiade for initial values of E, G, and MU. 

9. If E and >NU or G and NU are both blank, they will both be given the value 0.0. 

10. Implausible data on one or inore .'lATl cards will result in a warning -nessage. 
Implausible data is defined as any of E < 0.0, G < 0.0, 0.5 < NU < 0.0, or 

1 1 - i^ u' )G I ^ (except for cases covered by Remark 9). 
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Input Data Card .MATSl Material Stress Dependence 



Oescriotlon: Specifies table references and material properties which are stress-dependent for 
use in material nonlinearity applications. This card will be activated If MATl with the same MID 
Is being used In the material nonlinear solution sequence (66). 



Format and Example: 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


MATSl 


MID 


R1 


TYPE 


3 


YF 


HR 


LIMITl 






MATSl 


17 


28 


PLASTIC 




2 


3 


2.^A 




•^ABC 














LIMIT2 


















♦ABC 





















Field 

MID 

n 

. nPE 
3 



YF 



HR 



LIMIT! 

LIMIT2 



Contents 



Identification number of a MATl card (Integer > 0) 

Identification number of a TABLESl card (Integer > 0) 

Type of material nonlinearity (BCD) HLELAST (Nonlinear elastic) or PLASTIC 
(Plasticity theory) 

Wort hardening slope (slope of stress vs. plastic strain) In units of stress 
(Real). For elastic-perfectly plastic cases, B * 0.0. For more than a single 
slope In the plastic range, the stress-strain data must be supplied In a TABLESl 
card referenced In field 3. See Remarks 1 and 2. 

Yield function (Integer); one of the following values: (Default » 1) 

1 (von Mises) 

2 (Tresca) 

3 (Mohr-Coulomb) 

4 (Drucker-Prager) 

Hardening Rule (Integer); one of the following values: (Default ■ 1) 

1 (Isotropic) 

2 (Kinematic) 

3 (Combined Isotropic and kinematic hardening) 

Parameter used In the yield function specification (Real). See Remark A. 

Parameter used In the yield function specification for yield functions 3 and A 
(Real). See Remark 4. 
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Remar'cs: 1. If type ■ NLELAST, the stress-strain data given fn the TABLESl card will be use- 

determine the stress for a given value of strain. The values B, YF, HR, LIM 
and LIMIT2 will not be used In this case. 



2. If type • PLASTIC, either the table Identification, Rl, or the worit harder 
slope, 3, may be given but not both. If the table ID Is omitted, the i 
hardening slope, 8, specified in n-eio 5 Is defined as 



3 




or 




where E Is the elastic modulus and f’ Is the slope of the uniaxial stress-str 
curve Irr the plastic region. 

3. If Rl Is given, TABLESl entries of stress-strain data must conform to the follow 
rules: 



a. If type * NLELAST, the curve must pass through origin (x^»0., y^^O., for s 
value of 1). This option Is valid only for line elements. 

b. If type ■ PLASTIC, only data for the first quadrant need be supplied. Th 
are two notions available to specify the data. In case L (defined as F0RI 
In the TABLESl card) the first point Is at the origin, the second point Is 
the yield stress given In field 3 of the MATSl card. The slope of the 1 
joining the origin to the yield stress must be equal to the value of £ g1' 
on the MATl Bulk Data card. In case 2 (defined as F0RM»1 in the TABU 
card) the sloqes of the curve In the plastic range are supplied. See ' 
TABLESl card description for additional details. 

4. LIMITl and LIMIT2 are parameters used In the yield function definition as follows 



Yield Function 


LIMITl 


LIMIT2 


von Mises (1) 


Yield stress In tension, a 


Hot used 


or Tresca (2) 


y 




Mohr-Coulomb (3) 


Cohesion, c (In stress units) 


Angle of Internal 


or Orucker-Prager (4) 




friction ^ (In degree: 



S. See Section 2.14.1.2 In the MSC/NASTRAM Application Manual for a detail 
description of yield function and hardening rules. 
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Input Data Card TABLESl Material Property Table, Form 1 



Description: Defines a tabular function for stress-dependent material properties sucb as tne 

stress-strain curve and creep parameters 



Format and Example: 




Field 

ID 

I 

*1 



Contents 

Table Identification number (Integer > 0). 
Not used 

Tabular entries (Real) 



Remarks : 1. The must be In either ascending or descending order but not both. 

2. Jumps between two points (x^ » allowed, but not at the end points. 

3. At least two entries mist be present. 



4. Any x-y entry may be Ignored by placing the BCD string SKIP in either of the two 

fields used for that entry. 

5. The end of the Cable is Indicated by the existence of Che BCD string ENOT in either 

of the t'wo fields following the last entry. An error is detected if any 

continuation cards follow the card containing the end-of-table flag -NOT. 

6. TABLESl is used to Input a curve In the form of 



1 . Y^(X) 



where X is Input to the Cable and Y is returned. The cable look-up x = X, 

is performed using linear interpolation within the table and linear extrapolation 
outside the table using the last two end points at the appropriate table end. At 
Junp points the average yj(x) is used. 
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Inout Data Card PL3A02 Pressure Load on a Two-Dimensional Structural Element 



Descriotlon: Defines a uniform static pressure load applied to two-dimensional elements. Only 

QUAd^, ShtAk, TRIAD, or TRIAX6 elements may have a pressure load applied to them via this card. 



Format and Examole: 



1 


2 


3 


4 


5 


6 


7 


a 


9 


10 


PL3AD2 


SID 


P 


EID 


tID 


EID 


EID 


EID 


EID 


1 


PL3A02 


21 


-3.5 




4 


16 




2 




1 



Alternate Form; 



PL0AD2 


SID 


P 


EIDl 


"THRU" 


EID2 


^XlXj 






! PL3AD2 

i_ ! 


1 


30.4 


16 


THRU 


48 











Field 




Contents 






SID . 




Load set Identification number (Integer > 0) 






P 




Pressure value (Real) 






EID ) 
EIDl) 
EID2) 




Element Identification number (Integer > 0; EIDl < EID2) 






Remarks: 


1. 


EID must be 0 or blank for omitted entries. 








2. 


Load sets must be selected In the Case Control Deck (LDAD=SID) 
MSC/NASTRAN. 


to be 


used by 




3. 


At least one positive EID must be oresent on each PL3AD2 card. 








4. 


If the alternate form is used, all elements EIDl through EID2 
dimensional . 


must 


be t'wo- 




c 


The direction of the pressure is comouted according to the right-^’and rule using 
the grid point seauence specified on the element card. Refer to the PL3AD card. 




6. 


All elements referenced must exist. 








7. 


Continuation cards are not allowed. 







87 



Input Data Card NLPARH Parameters for Nonlinear Analysis Control 



Description : Defines a set of parameters for nonlinear analysis iteration strategy 



Format and Example; 



123456739 10 



NLPARM 


ID 


I.NC 


OT 


XMETH0D 


<STEP 


MAX I TER 


C0NV 


I.NT0UT 




NLPARM 


15 


5 




AUT0ON 






PW 




-NLl 




EPSU 


EPSP 


EPSW 


MAXDIV 


, MAXQN 


MAXES 


FSTRESS 






*NL1 


l.OE-3 


1.0-4 


1.0-5 


1 


20 


3 


.02 







Field 

ID 

INC 

or 

<METH0P 

KSTEP 

V63 MAX I TER 
C3NV 

INiaUT 

EPSU.EPSP, 

EPSW 

MAXDIV 

MAXQN 

MAXES 

FSTRESS 



Contents 

Identification Number (Integer > 0) 

Number of Increments (Integer > 0) 

Incremental time Interval for creep analysis (Real 2. (Default » 0,; no creep) 

Method for controlling tangent stiffness updates (3CD = AUT0, ITER .SEMI ,LSON,AUT0ON, 
or SEMIQN) (No default). See Remaric 5. 

Number of iterations before a .matrix update. Used only for ITER and '.SON 

methods. (Integer >^0 or blank) (Default » 1) 

Limit on total iterations for each load increment (Integer k 0) (Default 20) 

Flags to select convergence tests (BCD =• U, P, W, or any combination of the 
letters) (Default » W). See Remarks 3 and 4. 

Intermediate output flag (BCD =• YES or N0) (Default = YES). See Remark 5. 

Error tolerance limits for determining convergence (Real > 0.0) (Defaults » l.OE-2, 
l.OE-3, l.OE-4, respectively) 

Limit on diverging Iterations (Integer > 0) (Default = 1) 

Maximum number of Quasi -Newton correction vectors to be saved on the data base 
(Integer 2.0) (Default > 10) 

Maximum number of linear search operations per vector Iteration (Integer > 0) 
(Default » 3) 

Fraction of effective stress (T) used to determine the subincrement size in the 
material routines (0. < Real < 1.0) (Default » 0.15). See Remark 9. 



(Continued) 
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- ™- . 



^mm- 



P 



I 



NLPARM (Cont.) 



Rgnarics : 1 

2 

3 



. The NLPARM 3ulk Data card is requested by the Case Control card NLPARf^ = ID. Each 
solution subcase requires a loading condition and an NLPARM request. 

. In case of static analysis (OT » blank), INC is the numoer of equal subdivisions of 
the net load change defined for the SUBCASE. Aoolied loads, gravity loads, 
temperature sets, enforced displacements, etc., define the new loading condition. 
The differences from the previous case are divided by INC to define the incremental 
values. In case of creep analysis (OT > 0.), INC is the numoer of time step 
Increments. 

. The test flags (U.P.W) and the error limits (EPSU, etc.) define the convergence 
criteria (U » Olsolacement error test, ? = Load equilibrium error test, W = Work 
error test). If the internally calculated error fractions are less than the 

requested limits, the iteration stoos, the results are processed, and the program 
continues to the next load increment or subcase. All requested error criteria must 
be satisfied for a “convergent solution". 

If convergence does not occur for a particular load increment, the number of 
iterations is limited to MAXITER. If MAXITER is a positive number, the results 
will be processed and the program will continue with the next load increment. If 
MAXITER is a negative value, the solution will terminate if convergence is not 

achieved. 

. The basic nonlinear solution methods are the “QFGS" Quasi-Newton and .modified 
“Newton-Raphson" processes whereoy the out-of-balance nonlinear loads are measured 
and recycled to solve for improved displacements. The convergence of this 
iteration process is generally improved by either performing “line searches" or by 
frequent updates to the "tangent stiffness matrix" to account for current plastic 
and geometric effects. However, the matrix updates may be relatively costly and 

may be detrimental if the estimates are far from the correct answer. Thus, the 

user is given the following KMETH0O options to provide control over a large range 
of problem types. 

a. If the AUT8 option is selected, the program will automatically select the most 
efficient strategy based upon convergence rates. At each step the number of 
iterations required to converge is estimated. If this exceeds the '^AXITER 
limit, or if the CPU time exceeds the time to perform a stiffness matrix 
update, the module will perform a stiffness update. If divergence occurs for 
MAXOIV successive iterations, a stiffness update is also performed. 

b. If the SEMI option is selected, for each load increment the program will; (1) 

perform a single displacement iteration based upon the new load; (2) update 
the stiffness matrix for the estimated geometry and .material properties; and 
(3) resume the normal AUTO iteration method. 

c. If the ITER option is selected, the stiffness .matrix will be updated only 

after each <STE? iteration. The count is reset to zero for each new matrix or 
load step. 

d. If one of the AUT0QN, SEMIQN, or the LSQN options is selected the "Quasi - 

Newton" and “Line Search" procedures are used. Basically the line search (LS) 
will scale the displacement increment vector to minimize the error at each 
Iteration. If the line search provides a substantial improvement, the Quasi- 
Newton (QN) operations use the local tangent of the error to correct the 

tangent matrix solution for subsequent iterations. Otherwise the operations 
behave much like the AUTO, SEMI, and ITER options. 
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6. Output requests for ELFflRCE and STRESS made 1n Case Control will Pe processed if 

INTOUT ■ YES (tPe default value). If INT0UT ■ NO, tnese output requests are 
processed for only the last load factor In a subcase. Requests for OISP printout 
and structure plots are also processed during the nonlinear Iteration 1f PARAM, 
NLOISP, 1 Is provided In the 3ullc Data deck. 

7. If a diverging solution Is encountered for MAXOIV successful Iterations, the AUTO, 

AUTOQN, SEMI, and SEMIQN options will perform a tangent stiffness update. If, 

after the stiffness update, the solution Is still diverging, the program 
terminates. For ITER and LSQN options, diverging cases cause Immediate 
termination. 

3. The unit of the Incremental time interval under OT field nwst be consistent with 

the unit used on the CREEP card to define the creep characterl sti cs. Total creep 

time for the subcase may be found from OT multiplied by the value in the INC. 

9. The number of subincrements In the material routines (el astopl astic and creep) are 

determined such that the subincrement size Is approximately FSTRESS * a. 

FSTRESS is also used to establish a tolerance for el astopl astic material in the 
material routine, i.e.. 



6 in yield function < FSTRESS * yield stress 



If the limit is exceeded at the converging state, the program will exit with a 
fatal error message. Otherwise the stress state is adjusted to the current yield 
surface, resulting In 5 ■ 0. 
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Input Data Card $PCG Constraint Set Generator, Standard Form 



Description: Generates SPCl cards that constrain a subset of grid points In a grid point field 



Format and Examole: 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


SPCG 


SID 


FID 


C 


D1 


D2 










SPCG 


15 


3 


246 


0002 


0312 











Field 



Contents 



SID 



ID number of single-point constraint set ( Integer > 0) 



FID 



ID 'number of a GRIDG or GRIDU field 



C Component numbers (any combination of the digits 1-6 with no Imbedded blanks) In 

the displacement coordinate system (CD entry on GRIDG card) 

D1,D2 Location number pair of order-independent pivot points that define a subset of grid 

points In FID. ' See Part I, Section 2.5. 



1. See Part I, Section 5.0 for a discussion on the SPCG Input card. 



Remarks: 



2 . 



SPC entries on generated SPCl cards are made only for active grid points. See 
Part I, Section 10.1. 



Input Data Card 0L3AD Dynamic Load Comoinatlon (Superposition) 



Description : Defines a dynamic loading condition for frequency response or transient response 

proolems as a linear combination of load sets defined via RL3AD1 or RL3AD2 cards (for frequency 
response) or TL3AD1 or TUJAD2 cards (for transient response) 



Format and Examples: 



1 


2 


3 


4 


S 


1 


7 


t 


9 


10 


0L3AD 


SID 


S 


SI 


LI 


S2 


L2 


S3 


L3 




DL3AD 


17 


1.0 


2.0 


6 


-2.0 


7 


2.0 


3 


-A 






S4 


L4 




-etc,- 












‘A 


-2.0 


9 

















-etc.- 



Field 

SID 

S 

Si 

Li 



Contents 

Load set identification number (Integer > D) 

Scale factor (Real) 

Scale Factors (Real ) 

Load set identfication numbers defined via card types enumerated above 
(Integer > D) 



Remarks: 1. 



The load vector being defined by this card is given by 



iP} “ S I S.lP-i . 



2. The Li must be unique. 

3. SID must be unique from all Li. 

4. Nonlinear transient loads may not be Included; they are selected separately in the 
Case Control Deck. 

5. Linear load sets jrust be selected in the Case Control Deck (DL3AD = SI0) to be used 
by MSC/NASTRAN. 

6. A DL3AD card may not reference a set identification number defined by another 0L3AD 
card. 

7. TL3AD1 and TL3AD2 loads may be combined only through the use of the 0L9AD card. 

8. RL3AD1 and RL(JAD2 loads may be combined only through the use of the DL9AD card. 

9. SID must be unique for all TUJADl, TL3AD2, RL3AD1 and RL3AD2 cards. 
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! 
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Input Data Card TL3AD1 Transient Response Dynamic Load, Form 1 



Description : Defines a time dependent dynamic load or enforced motion of the form 



lP(t)| • (A f(t - t)( 



for use in a transient response problem 



Format and Examole; 



1 


2 


3 


4 


S 


6 


7 


3 


9 


10 


TU7ADI 


SID 


L 




TYPE 


TID 










TL3A01 


5 


7 


9 


0 


13 











Field 

SIO 

L 

H 

TYPE 

TID 



Contents 

Set identification number (Integer > 0) 

Identification number of DAREA card set or a static load card set or a thermal load 
set (in heat transfer analysis) which defines A (Integer > 0) 

Identification number of DELAY card set which defines t (Integer 2, 0) 

Defines the nature of the dynamic excitation (Integer 0, 1, 2, 3 or blank), (See 
Remark 2.) 

Identification number of TABLEDi card which gives F(t - tJ (Integer > 0) 



Remarks : 1. If M is zero, t will be zero. 



2. The nature of the dynamic excitation is defined In Meld 5 in accordance with the 
following table. 



Integer 


Excitation Function 


0 or blank 


Force or Moment 


1 


Enforced Displacement 


2 


Enforced Velocity 


3 


Enforced Acceleration 



See Section 2.7 of the MSC/NASTRAN Application Manual regarding the use of the 
enforced motion options. Note that “large masses" must be used for enforced 
motion. For heat transfer problems. Field 5 must be blank. 

(Continued) 
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TL3A01 (Cont.) 



3. Dynamic load sets mist be selected In the Case Control Oec)( (0L3A0*$I0) to be used 
by MSC/NAS7RAN. 

4. TL0AOI loads may be comoined witn TL3A02 loads only by soeci f 1 cati on on a 0L3A0 
card. That Is, the SIO on a TL3A01 card may not be the same as that on a TL3AD2 
card. 

5. SIO must be unique for all TL3A01, TL0AO2, RL3A01 and aL9A02 cards. 

6. Field 3 may reference sets containing QHBOY, QBOYl, QB0Y2, QVECT, and QV0L cards 
wnen using the heat transfer option. 

7. If the heat transfer option is used, the referenced QVECT data card may also 
contain references to functions of time, and therefore A may be a function of time. 

3. Fourier analysis will be used if this is selected in an aeroelastic response 

problem. 
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Inout Data Card LSEO Static Load Set Definition 



Description; Defines a seouence of static load sets to be applied to the structural -nodel . The 
toad sets may be referenced by dynamic load cards. 



Format and Example; 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


LSEO 


SID 


OAREA 


LID 


TID 












LSEO 


100 


200 


1000 


1 1001 













Field Contents 

SID Set Identification of the set of LSEQ cards (Integer > 01 



DAREA 



The DAREA set identification assigned to this static load vector (Integer > 0) 



LID 

TID 



Load set Identification number of a set of static load cards (Any card that may be 
referenced by the L3AD Case Control card) (Integer > 0 or blank) 

Set Identification of a thermal load set (Any card that may be referenced by the 
TEMP(L3AD) Case Control card) (Integer > 0 or blank) 



Remarks: 1. The above cards will not be used unless selected in the Case Control Deck with a 

L3ADSET card. 

2. This card is available only in suoerelements and in Solution Sequences 25, 27, 30 
and 31. 



3. The number of static load vectors created for each suoerelement is the numbe'- of 
unique OAREA IDs on all LSEQ cards in the bulk data. 

A. The OAREA Identification assigned to the static load vectors may be rofpronced by 
RL3A01, RL3AD2, TL3AD1 and TL3A02 cards. 

5. Element data recovery for thermal loads is not currently implemented in dynamics. 

6. OAREA set identification numbers should be unique with respect to all static load 
set identification numbers. 



0 
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Input Data Card TABLEOI Dynamic Load Tabular Function, Form 1 



Description : Defines a tabular function for use In generating frequency-dependent and tlme- 

dependent dynamic loads 



F ormat and Example: 



12 3 A 5 6 7 3 9 10 




- etc.- 



Fleld Contents 

ID Table Identification number (Integer > 0) 



*1 *-yi 



Tabular entries (Real) 



Rema r>cs : 1. The must be in either ascending or descending order but not both, 

2. Jumps bet'^een two points (x^ * + allowed, but not at the end points. 

3. At least two entries mist be present. 



4. Any x-y entry .may be Ignored by placing the 3C0 string ''SKI?” in either of the two 
fields used for that entry. 

5. The end of the table Is Indicated by the existence of the 3CD string "iNDT" in 
either of the two fields following the last entry. An error is detected if any 
continuation cards follow the card containing the end-of-table flag "ENDT". 

6. Each TABLEDI mnemonic Infers the use of a specific algorithm. For TABLEDl type 

tables, this algorithm Is 



Y > yt(x] 



where X Is Input to the table and Y is returned. The table look-up yj( <j , x = X, 
is performed using linear interpolation within the table and linear extrapolation 
outside the table using the last two end points at the appropriate table end. At 
juitp points the average y-p(x) is used. There are no error returns from this table 
look-up procedure. ‘ .. 

7. Linear extrapolation Is not used for Fourier Transform methods. The function Is 
zero outside the range. 
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Input O^ta Car'd T ZlZ?*iL Step an^j Nonlinear Solution Ccntrolt 

Description : Provides pararnetric controls ano Cata 'or Nonlinear 'transient ^nal/tlt 

romat and i;fafTOle: 



1 


2 


3 


4 


5 


5 


7 


z 


9 


10 


T3TEPNL 


ID 


•®T 


3T 




.“PTrOO 


f.DTE? 


.•'A/. IT'?. 


33'iV 1 


tstepnl 


250 


40 


.3015 1 


i 5 i 


1 4'JT2 




5 


3 i 






£?SU 


iPZ? 


£?r- 


“Ji/OIY 


!^A/CN 1 




} ^DT^ESS 






^NLl 




.01 






23 1 


iS 


.02 


i 



Field 

ID 

not 

OT 

N9 

‘-ETH30 

<3TEP 

C3NV 

t?SU,£?S?.i 
.^AXDIV 
■‘-.AXON 
?^AXLS 
' 5TPE3S 



Contents 

Identl ficatlon nurr.oer (Inte^ger > 0) 

Numcer of time steps of /alue 01(1) (Intecer > 2) 

Tine Incrernent (Peal > 0*0) 

£<1p factor for o^jtput 'every Ngtn step ^lll te saved 'or :;jtout) ^Intec?'“ > 0) 

‘'ethod for controlling Iterations and tangent matr*/. updates or "ST£?) 

(No default ) 

Nurncer of steps tefore a matrix update (Irteger > 1) (Default * C) 

Limit on total Iterations 'or eacn tl'^e Increr’.ert 'Integer / C) 'Sefauls * 2) 

Flags to select convergence tests (sCD * J, or any cor^'natlon of tre 

letters or olanic) 

PTrt zrror tolerance Units 'or detertlning convergence '^eai y D.D) 

(Defaults » l,S£-2, l.CE^S, l.CE-A, respect*' vely ) 

Limit on tne nijrcer of d1/erg1ng Iterations at eacn tire and on ol/ero'^c t^me 
steps (Integer > 1) (Default » 2) 

J^axinun nuraer of Duasi -Nekton correction /ectoro to oe saved on tre data case 
(Inrecer >_ 0) p?fi‘jtz » 20) 



?*2xi3un r.unze'" of lire Ssiron i^erji’cr^s *or 
(Oe^i'jlt * 2) 



'.;:n x-.'cr '.i-.ic*' > 0) 



''■jctlon of *f^»ct 1 vo stress cj js?c *3 tri* '331 '"2% *■ 

T.ater*;al 'ouclnes (0. < 2»?1 < 1.3) » 3.1;). ;e« ’,0. 



{Zor.Zir-jiO ) 
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NASTRAN DATA DECK 



TSTEPNL (Cont.) 



Refflartcs : 1. The TSTEPNL Sulk Data card Is selected by the Case Control card TSTEPNL • ID. Each 

subcase (residual suoerelement solutions only) requires a TSTEPNL card and either 
aooHed loads via TL3A01 data or Initial values from a previous subcase. Multiple 
subcases are assumed to occur sequentially 1n time with the Initial values of time 
and displacement conditions of each subcase defined by t.he end conditions of the 
previous subcase. 

2. OT, SOT, NO are the time Increments, At, the total number of steps N^, and the out- 
put Interval, ng, respectively. The value of time, T^, at each Increment n will be 



tn ■ tg -*■ ndt n » 0, l,...,N, 



where tg Is the Initial value of time which corresponds to the last step from the 
previous subcase. 

3. For printing and plotting the solution, data will be output at steps n » 0, ng, 

2ng, .N. The Case Control card OTIME may also be used to control the output 

pol nts. 

A. The basic solution method Is the "Newmarlc Beta" transient numerical method along 

with the "Newton-Raphson" method for controlling nonlinear effects. Equilibrium 
solutions with respect to Inertial, damping, elastic, and nonlinear loads are 
obtained at each time step much the same as In the statics case. 

The KMETHOO options provide the user with three optional nonlinear iteration 

methods. These are: 

AUTJ provides for automatic stiffness matrix updates to Improve poor 
convergence. The KSTEP value Is Ignored. 

TSTEP provides for a matrix update at every KSTEP^^ increment of time 
(typical ly a large number). 

5. The stiffness matrix Is always updated for a new subcase or restart. Only one 

update is allowed for each time step and the Iteration process returns to the 

previous time step after an update. 

6. MAXITER defines the number of static nonlinear load Iterations (correcti ons ) to be 
performed at each time step. At least two Iterations are necessary to determine 
the convergence of the displacements (see below). 

7. The test flags (U,P,W) and the error limits (EPSL), etc.) define the convergence 

criteria (U » Displacement error test, P » Load equilibrium error test, W » Uor\ 
error test). If the internally calculated error fractions are less than the 

requested limits, the iteration stops, the results are processed, and the program 
continues to the next time step or subcase. 
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TSTEPNL (Cont.) 



9. If 10 tests are requested or no convergence occurs, tne numoer of Iterations is 

limited to MAXITER. The results *<111 be processed and the program <*111 continue 
with the next time step. 

9. MAXQN provides the control over diverging solutions. At each time step, the 

solution Is allowed to diverge until, at the MAXOIV Iteration, a matrix update Is 

requested. If the average convergence parameter"^ Indicates divergence for ;<AXOIV 

successive time steps, the transient solution stops and exits for data recovery. 

10. The number of subincrements In the material routines (elastoplastic and creep) are 

determined such that the subincrement size Is approximately FSTRE3S *7. 

FSTRESS Is also used to establish a tolerance for elastoplastic material In the 
material routine, i.e.. 



5 in yield function < FSTRESS • yield stress 



If the limit Is exceeded at the converging state, the program will EXIT with a 
fatal error message. Otherwise the stress state Is adjusted to the current yield 
surface, resulting In 5 =■ 0. 
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LINEAR CEEORHATIOK WITH OYAAKIC LOADS OECEHaER 1. 1983 HSC/NASTRAN b/ 1/83 PAG 

ANSIEN1 RESPONSE EV SNCCK kAVE LOAC. 
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INPLI BULK OAT A CARD COUNT 



kCNLlKEAR DEFORMATION hllH OTNAMIC LOADS DcCEMBER 1> 1963 MSC/NASTKAN Ml 1/B3 PAG 

TRAKSIEM RESPONSE EV ShCCK kAVE LOAC. 
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transiem response ev shcck lave loac. 
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TRANSIENI RESPONSE EY SHCCK SAVE LOAD. 
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